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A METHOD FOR DETERMINING SUPPLY QUANTITY FOR 
THE CASE OF POISSON DISTRIBUTION OF DEMAND* 


M. R. Mickey! 


General Analysis Corporation 
Los Angeles, California 





The problem treated is that of determining the number of units of an 
item to be supplied against a future demand, on the basis of the past 
record of demand for that item. A method is described for the case of 
Poisson distribution of demand and known relation between the expected 
demand for the experience and the future periods. 











INTRODUCTION 

Inventory problems are sufficiently complex that flexible, relatively simple, non- 
optimal solutions may be of practical as well as theoretical interest. The problem of present 
concern is that of determining the number of units of an item to be supplied against a future, 
unknown demand. The emphasis here is on that aspect of the problem that arises when the 
distribution of demand is not known but must be estimated on the basis of experience., The 
characteristic of the solution that is of interest is that the expected value of the performance 
characteristic will not exceed a preassigned bound, whatever the values of the unknown 
parameters of the demand distribution may be. 

The discussion is confined to a single solution for a particular case. The case is 
interesting in itself, and the treatment is sufficiently suggestive to indicate extensions. In 
particular, the following problem is considered. Consider the routine supply of a single spare 
part. Suppose that X units are requested over a given period of experience. It is desired to 
supply a quantity u(X) units against the demand of a specified future period. The problem is 
to choose the function u(x) such that the probability of failure to meet the demand does not 
exceed a given bound, say @,. 

The requirement is ordinarily interpreted as applying it to the conditional probability 


(1) P{Y>u(X)|X=x)}=a 


where Y is the random variable denoting future demands and Xo is the particular value 
observed for the experience period. This requirement is overly restrictive and in general 
precludes the existence of a solution. A second interpretation, that adopted here, is simply 
that the inequality 





*Manuscript received September 8, 1958. 
Prepared while the author was at the RAND Corporation as P-1376, 15 May 1958. 
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(2) P{y>uwx)}<a 


be satisfied for all values of the parameters of the demand distribution. This interpretation 
is consistent with the relative frequency approach to probability. In a series of N trials the 
expected number of times that demand exceeds supply is not greater than Na. Note that this 
is the case for consecutive trials, for which what is future demand for one period becomes 
part of the past experience for the next. 

Explicit solution of the problem of satisfying inequality (2) will depend upon the assumed 
functional form of the demand distribution. The case presented here is that in which X and Y 
are (independently) distributed according to Poisson distributions with means k,A and k »; ky 
and are assumed known and A is unknown. One knows only the relation between expected 
demand in the future and the experience periods. A solution to this case is obtained by observ- 
ing the following equalities: 


P{ y>u(x)} = P{x+¥Y>X+u(x) } 
(3) 
E{P{z>x+ux)|x+¥ = z}} 


and by noting that the conditional distribution of X given that Z = z is the binomial distribution 
with parameters z and p, [1, 2],! where 


(4) — a 
~ + 


Let z(x) be the smallest integer z such that 


x 
(5) >: (7) pla-n*) <e, 


j=0 


and define, for integral values of x, 
(6) ug(x) = 2(x) - x - 1. 


Let x(z) be the function inverse to z(x) for z = z(0). Then, since z(x) is an increasing 
function of x, 


(7) X + up(X) = 2(X) - 1< Z 
if and only if both 


(8) X = x(Z) 





lThese references were kindly supplied by D. S. Stoller, who has independently noted the idea 
developed in this paper. 
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and 
(9) Z = 2(0) 


are satisfied. Since the definition of z(x), Equation (5), is equivalent to 
(10) P{X =x(Z) | Z}<a, Z =2(0), 
it follows from (7) - (10) that 
(11) P{x+ug(x) <z|z}<a. 
This combined with (3) gives the desired result 
(12) P fy > ug(x)| sa, 
and hence the function Ug(x) defined by (6) is a solution to the problem. 
For small value of x, values of u(x) can be obtained from various tables. For example, 
for given p, tables of the binomial distribution may be used to obtain the smallest value of z 
satisfying inequality (5). Alternatively, the sum of inequality (5) may be evaluated by use of 


tables of the incomplete beta function. In this form, u(x) is obtained as the smallest integer u 
such that the inequality 


1 
(13) il —Mlu+x+2) 4% (1-4)" at < a 
p T(u+l) I(x+1) 


is satisfied. If, in this integral, the transformation 


(14) " (u+1) t 
(x+1) (1-t) 

is applied, the integral of (13) is the upper tail probability of the F distribution, for parameter 

values 2x+2 degrees of freedom for the numerator mean square, and 2u+2 degrees of freedom 

for the denominator mean square. Therefore (13) is equivalent to the inequality 


(15) u + 1 = (x+1) — F , (2x+2, 2u+2) 


where F, (v1, Vo) is the 100a percentage point of the F distribution with Vy and Vo degrees 
of freedom. It is of interest to note that since x/k, is an estimate of the unknown A, 


(x+1) (1-p) (x+1)k ‘ 
(16) - hed 


p ky y 





is an estimate of a, = kA. In these terms inequality (15) may be expressed as: 
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(17) u+12 y Fy (2x+2, 2u+2) , 
(18) (u+1) Fy_. (2u+2, 2x+2) = ly , 

The form of the inequality given by (18) is interesting in that it leads directly to an 
important property of the solution. This is the property that as experience increases, the 


solution converges to that for the case of known A. Formally, let x and ky tend to infinity 
such that y converges to Ay Then inequality (18) becomes, in the limit: 


1 
(19) (us) Fy_q (2042, @) => x? (2u+2) = dy 


where - (2u+2) is the 100(1-@) percentage point of the chi square distribution with 2u+2 
degrees of freedom. The solution for known A y is given as the smallest integer u such that 


(20) 


Since the sum of (20) can be evaluated as a chi square integral, 


2u+2 _ 





, cf V 
set (3) 
Aye -{ 7% 
j! 


u+1 0 


inequality (20) is equivalent to the requirement that 2A_ be not greater than the percentage 
point of the chi square distribution, (2u+2). This is the requirement of (19), and provides 
the central point of the demonstration. 

For large values of x and u, the solution can be approximated by applying the normal 
distribution approximation to the binomial sum of (5). In this case z is the smallest integer 
not less than z' where 


1 
(22) x +; = pz'-t, Vz' p(1-p), 


and t, is a percentage point of the normal distribution. Equation (22) may be solved for 
u'=z'-x-1as 





t 
" l-p ,.2 ay) | 2 | 
23 ee - Bis =e ” a 
(23) +5, “a tz ~~ *<; (1-p) t7,-2 
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where, as before, 


os 


k 
A. = La yon = —Y (x+1) ‘ 
p k 


y 
x 


If the second term in the radical is neglected, the approximation becomes the integer nearest 
u" where 


2 
. ky (2-1) ty ¢ 
(24) ut =f} +7 \2%@ / 4 tyr ‘ 
bf 2k, Np b | 


An interesting point with respect to (23) and (24) is that these exhibit the "additional 
stock" needed to compensate for the uncertainty in the mean demand. If (24) is expressed in 
the form 


2 - 
(25) u" = lay +ty Viy |+ ae (- 1) to yiy 


the first term on the right hand side represents the stock that would be required if the mean 
demand were known and equal to d , and the second term represents the additional stock 
required to compensate for the uncertainty. If the experience base is large relative to the 
prediction period, that is, p is close to unity, then but a small percentage increase in stock 
is required. 

This same point can be examined for small value of u by use of (17). Table 1 gives a 
fragmentary table of values of \_ for a = 0.05, and for some values of x and u. The stockage 
function u may be regarded as a function of }_ and x, u=u(4i_, x). Fora given x, u isa 
step function of 1 the jumps occur at values of A_ for which equality holds in (17). It suf- 
fices to tabulate these values of A, and this is the basis for Table 1. For example, for x= 0, 
i.e., no demand during the experience period, 


u=- 0, Ay = 0.0526, 
u=1, 0.0526 < , = 0.288, etc. 


The interesting point brought out by Table 1 is that for small values of u, the solution is very 
nearly that given by assuming X,, to be the true demand. 

Although the preceding development depends upon a special property of the Poisson 
distribution, it is clear that the underlying idea is not restricted to the case of Poisson demand 
distributions. Also, the performance criteria need not be limited to that of the probability of 
meeting the demand. For example, the number of shortages could be taken as the character- 
istic of interest, and the number of items to be supplied could be determined so that the 
expected value of this quantity would not exceed a prescribed bound, regardless of the values 
of the parameters of the demand distribution. 
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TABLE 1 
Values of y = Ay (u,x) for which jumps of u(A_; x) occur, | (16) 
for 95-percent assurance of meeting demand; u(A y’ xy=u 

for A (u-1, x) < ry = A,(u, x). 






































u 
0 1 2 3 4 5 
Tab 
= befc 
0 0.0526 | 0.288 | 0.584 | 0.897 | 1.220 | 1.542 _ 
5 0.0515 | 0.338 | 0.750 | 1.220 | 1.718 | 2.230 
10 0.0514 | 0.346 | 0.778 | 1.276 | 1.816 | 2.377 a to 
100 0.0513 | 0.354 | 0.813 | 1.356 | 1.953 | 2.586 
oc 0.0513 | 0.355 | 0.817 | 1.365 | 1.968 | 2.608 
Solutions along the lines of this paper are of interest in that they provide a flexible 
procedure for compensating for the uncertainty in values of parameters of the distribution of 
demand. In addition to possible practical value, these nonoptimal solutions may be of interest 
in evaluating the degree to which this uncertainty is important. 
ILLUSTRATIVE EXAMPLES 
Example 1 
Suppose that in a six months' experience period, 4 wings of aircraft had zero demand Thu 
for a given part. How many units of the part should be supplied for 6 wings for a three month the 
period for 95 percent protection level against shortage? = 
Let A be the mean demand per wing month. It is assumed that the experienced demand | this 
is an observation on a random variable distributed as the Poisson distribution with mean 24), m= 6ine 
and that the demand for the future period is a Poisson distributed random variable with mean 
18; k, = 24, k, = 18. mB (19) 
Consider first the direct method of solution based upon (4) and (5) ; 
is s 
i ae 
 —_— =" G OExz 
k, + ky 42 7 
sta: 
Since @ = 0.05, z is determined as the integer satisfying oa 
Zz z-1 : ass 
(l-p)*<a@ <(1-p) a 
z z-1 : 
(3) = 0.05 < (=) : ry | 
35- 
; at t 
The inequality is satisfied for z = 4, and hence u= z-x-1=4-0-1= 3, so that three 
units should be supplied for 95 percent protection against shortage. Note that if p= 1-a, no tha 


units are to be supplied. 
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The example may also be solved by use of Table 1. For this purpose compute hy from 
(16) as 


k 
5 = (met) = 2 ost) = 0:78. 
yk 24 


Table 1 is entered with x= 0, and since 0.584 < 0.75 = 0.897, u(r = 0.75; x = 0) = 3, so that, as 
before, 3 units are to be supplied. 


Example 2 
Suppose that in the above example, the experience had been 2 units per wing-month, i.e., 
a total of 48 units. Consider the approximation given by (24), where t, = 1.645 for @ = 0.05. 


k 
en = 29) (9) « 36.75 


2 (24) 


~» 
ul 


2 
u = (36,75) + (28) [(1.645)"-1] | 1.645 eae 


(2) (24) ‘ [4 
7 


(36.75) + (0.64) + (13.19) = 50.58 . 





i} 


Thus, 51 units should be supplied. This is the result given by linear interpolation in tables of 
| the percentage points of the F distribution. Numerically this interpolation is not sufficiently 
accurate to satisfactorily distinguish between u = 50 and u = 51. It is of interest to compare 





d this value with that obtained if Ay were known equal to the estimated value, 36.75. The 
inequality 
b 
(19) x? | (Qu) << 2a. <X2 . (2u42) 
l-a y~ * l-a 


is satisfied, in this case, for u = 47. 


Example 3 
a Suppose that the parts required in Example 2 are to be delivered in monthly intervals 
starting with the beginning of three months period. How should the delivery quantity be 
scheduled so as to maintain a 95 percent protection against shortage for each month? It is 
assumed that if shortage occurs during the first or second month it will be filled, if possible, 
from the next delivery. 
: For the first month, ky = 6, AL = 12.25, u= 19. For the first two months, ky = <2, 
A, = 24.50, u = 35. Hence, the appropriate delivery schedule is: an initial supply of 19 units; 
35-19 = 16 units to be supplied at the end of the first month; and 51-35 = 16 units to be delivered 
at the end of the second month. 
A more involved form of the problem is that of determining the delivery schedule such 
that there is 95 percent protection, say, against shortage throughout the three month period. 
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This form of the problem is beyond the development of this paper, although the same approach 
can be applied. 
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PROLEGOMENA TO OPTIMUM DISCRETE SEARCH PROCEDURES* 


Nelson M. Blachmant? 
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An object may appear in any one of n locations and, having appeared, 
will remain there. The probability of appearance in the ith location is 
pj, and the time of appearance is distributed uniformly over a long 
interval. A look in the ith location takes a time t; and it detects the 
object with probability P; if it is there. A search procedure, i.e.,a 
sequence of locations in which to look, is desired which minimizes the 
expected delay between the appearance and the detection of the object. 


It is shown that the optimum search pattern ideally involves regularly 
spaced looks in each location, the frequency of looks in the ith location 
being proportional to 


If P; is a function of t; and the latter is chosen so as to minimize the 
expected delay, then, when n is large, 


: 


dt; t; 


Ifthe probabilities p; are unknown, it is desirable to know what set of p; 


maximizes the expected delay. The worst p; is proportional to 
ti fe - 2)" and the corresponding optimum frequency of looks in the ith 
i 


a : a; | 
location is proportional to (5. ~ 3); 


In the Appendix consideration is given to the case where the various 
looks in the ith location have different durations, but the work on this 
problem is incomplete. 











*Manuscript received August 5, 1958. 
TWork done under U. S. Army Signal Corps contract. 
fon leave from Sylvania Electronic Defense Laboratory. 
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LIST OF SYMBOLS 
n the number of different locations in which the object may appear. 


Pj the probability that the location in which the object will appear is the ith. 
ey the duration of the search. 
t; the duration of a look in the ith location. 


P, the probability that a look in the ith location will detect the object if it is there 
from the beginning of the look. The look is assumed not to detect the object if it 
appears after the beginning of the look. 


Q the probability that a look in the ith location will not detect the object if it is there 
from the beginning of the look. Q; = 1- P; ‘ 

Jj the total number of looks in the ith location during the course of the search. 

Thy the interval between the beginning of the (j-1)st look in the ith location and the 
beginning of the jth. The zeroth look is defined to occur at the beginning of the 
search. 


t the expected delay between the appearance and the discovery of the object. 


INTRODUCTION 

Koopman [1, 2,3] has treated the problem of optimum allocation of a given amount of 
search effort to the discovery of an object known to be somewhere in a given region. The 
present paper treats a somewhat different problem, in which the object is not present at the 
beginning of the search but has a distribution of arrival times, and the aim is not to maximize 
the probability of detection but to minimize the expected delay between arrival and detection. 
Thus, the sequence in which to look in different places becomes important, and it is necessary 
to return, from time to time, to places already searched. Our problem may be stated as 
follows. 


An object may appear in any one of n locations and will thereafter remain there, the 
probability of the ith location being Pj, with 


n 
2m =}. 


i=1 


The time of appearance of the object is distributed uniformly over a very long interval of 
length T. A look in the ith location takes a time ti and, if the object is there, the look detects 
it with probability P,. In what order should the various locations be scanned during the time 
T to minimize the expected delay between the appearance of the object and its detection? We 
assume that if the object appears after the beginning of a look, it will not be discovered by that 
look. 

Any search pattern is completely characterized by the intervals Ti between the 
beginning of the (j-1)st look in the ith location and the beginning of the jth. Ti is defined as 
the interval between the start of the search and the first look in the ith location, it being 
assumed that the object will not appear before the start of the search. For a given search 
procedure, the expected delay between the appearance and the discovery of the object is 
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J, J-i 
Til 1 a 
‘ ij 
(1) t= 28 oe vr Tiga) | (> 
= }= = 


where J; is the total number of looks in the ith location and Q =1- P; is the probability of 
non-detection on each look. (If the object should go completely undetected, the associated delay 
would be the time from its appearance until the end of the search. However, when T is large, 
this event has a negligible probability.) 

To see where this expression came from, we note that (1) is given by the average, over 
all locations, of the expression in braces, which is the expected delay if the object appears in 
the ith location. The first term within the braces, ti, accounts for the duration of the look 
that discovers the target. The second term within the braces describes the expected delay 
between appearance and the start of the successful look. It, in turn, is the average, of the 
expression within the brackets, over all of the intervals between looks, during which the object 


T. 
can appear; the factor = is the probability of its appearance in the jth interval. The first 


1 
term within the brackets, 2 Tij? is the expected delay between the appearance of the object 


and the beginning of the first look in its location, and the second term accounts for the addi- 
tional delay due to the possibility of non-detection. The probability that the delay between 
appearance and detection will include Tigjat is Q;, the probability of mre on the 
first look; the probability that this delay will also include Ti(j+2) is Qa , the probability of 
non-detection on both of the first two looks; and so forth. 


DETERMINATION OF OPTIMUM Ti 
We now want to determine the characteristics of the search procedure that minimizes 
(1). Thus, we want to choose positive quantities Tj that minimize (1) subject to the conditions 


Ji 
(2) i.%** so 
j=1 


which simply mean that the last look in the ith location occurs about T after the start of the 
search process. In solving this problem we shall regard (2) as a strict equality, although it is 
only approximate; the resulting error vanishes as T becomes infinite. First we shall minimize 
(1) subject to (2), regarding the J, as fixed; then we shall determine the optimum J; subject to 
the condition that 


n 


(3) > 3% =T7, 


i=1 


i.e., T is the total time available for looks. Later we shall consider the additional constraint 
that looks in different locations must not overlap, i.e., we can look in only one place at a time. 
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Thus, from the right-hand side of (1) we subtract L; (a Lagrange multiplier) times the 
left-hand side of (2), summed over i, and we set equal to zero the partial derivative of the 
resulting expression with respect to every Ti? getting 


J; 
1 
‘ TL, 
(4) >, a) - Tt, = — (j = 1,2,...,355 ee ee 
k=1 Pj 


after some simplification. Solving these simultaneous linear equations for Tik? we find (by 
matrix inversion (see (A-2) and (A-3) in the Appendix)) 


1-Q, TL 
?., = STM (2 <j <J,-1) 


N 14Q Bj 





(5) 
Tk, 
: — (j =1, J;) ? 


1+Q Pj 





i.e., except for the last look, all looks in the ith location should be uniformly spaced. The 
Lagrange multipliers L; must be determined by substitution of (5) into (2). Thus, on the 
assumption that the duration of the search T is great enough that all J j are large, we find 





(6a) Ty = = (2 =j <J,-1) : 
Jj 
. 
(6b) = acter O08 1,4), 
(1 - QJ; : 





Because J i is large, the fact that (6b) is different from (6a) may be overlooked in 
determining the expected delay between the appearance and the detection of the object. Thus, 


putting Tij => in (1), we find, on carrying out the inner summations and dropping terms of the 
i 


T 
order of = (or changing the upper limit on the summation over k in (1) to infinity), that 


Jj 


2 T/1 1 
(7) t =)‘ », 5 (&- 2) 


i=1 


Equation (7) is easily seen to be correct without reference to (1); the first term within the 
brackets represents the duration of the look that discovers the object. In the second term, 


> represents the expected number of looks in the ith location while the object is there before 
i 
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itis discovered. From this J i must be subtracted because the object does not appear t 
P. 2 a, 
i i 


before the first of these looks but only > before, on the average. 
i 


OPTIMUM Ji AND MINIMUM t 
To determine the optimum J;, we must minimize (7) subject to (3). Hence, to the right- 


hand side of (7) we add L (a Lagrange multiplier) times the left-hand side of (3), and we set 
equal to zero the partial derivative of the resulting expression with respect to every J ip getting 


(8) J; x “i(4-1). 


The Lagrange multiplier is equivalent to the constant of proportionality in (8), which can be 
determined by substitution into (3). Thus, 


D, 
T —(>-3] 
& \2, 2 


(9) J, = a. ' 
n 
pan Pte (~~ ' =) 
i'=1 “+ * 


whence the optimum interval between successive looks in the ith location is 











&. 
- 





(10) T.. 


Thus, when the search is long, the optimum search pattern, described by the Tj,» becomes 
independent of the search duration T. This search pattern is optimum not only for finding a 
single object whose arrival time is distributed uniformly over T but also for finding many 
objects whose arrival times form a Poisson process (since the arrival time of each would be 
distributed uniformly)—unless they should arrive so frequently that the question of detecting 
two simultaneously would arise. Substituting (9) into (7), we find the minimum expected delay 
between the appearance and the discovery of an object, 


2 


(11) oe be a wa 5) 


i=1 


which is also independent of T. 
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Overlapping Looks 

In general, it is not possible to arrange a search pattern that satisfies (10) for all 
locations because we have not taken into account the condition that looks in different locations 
must not overlap. Thus, we can conclude only that a search pattern which approximately 
satisfies this equation is, at least, approximately optimum. However, it can be shown [4] that 
a non-overlapping search pattern based on (10) is asymptotically optimum, i.e., that, as the 
number n of different locations grows infinite and the probability of each approaches zero, it 
becomes possible to set up a search pattern with intervals given approximately by (10) for 
which the ratio of the expected delay to (11) approaches unity. (""Approximately" here means 
“with a ratio approaching unity.") 

If one confines his attention to periodic search patterns, i.e., to patterns which repeat 
after some arbitrary time interval, it is possible to treat a more general search problem, in 
which one seeks to maximize the expected gain instead of minimizing the expected delay, with 
the gain resulting from detection in each location given by an arbitrary non-increasing function 
of the delay between the appearance of an object there and its detection. A cost of a look in any 
location can also be included. 


Case of Arbitrary ti 


We may now suppose that tis the time required to look in the ith location, is arbitrary 
and that P;, the detection probability, is a function of ti. Thus, setting the derivative of (11) 
with respect to every ti equal to zero, we find that if we can ignore the first term of (11) (i.e., 
if ti is very much smaller than the expected delay), the optimum ti is that for which 


(12) ans 


This problem is considered further in the Appendix. 


Case of Unknown Pj 


If the a priori probabilities of the various locations p; are unknown, it is of interest to 
find the worst Pi» i.e., the values that maximize (11), so as to obtain an upper bound on the 
expected delay and so as to determine how to search in this case in order that the expected 
delay may be no larger than necessary. If, again, the first term of (11) can be ignored, we 


n 
find the worst set of Pj by maximizing (11) subject to the condition - P; = i, Viz, 
i=l 


(13) P, «t; (—-2). 
P, 2 


From (10), then, we have 
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and the expected delay is 


n 
(15) ‘<2. me 5) te 
i=] 
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APPENDIX A 

The Case of Arbitrary ti 

We may consider a more general case of arbitrary ti than that discussed above. 
on 


Instead of asstiming that all looks in the ith location have the same duration ti + a the same 
po of success P; =l1- Q; if the object is there, we nee ene that the ie look in the 

» location has pa, ti and qnsedinn probability Pi = 1- Q. In this case the expected 
he between the een and the detection of the object is no longer (1) but is 


ny 


y J; J;-j k-1 
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To minimize (A-1) subject to (2), we proceed as before. Instead of (4) we now get a set of 
simultaneous linear equations in the Tix in which the matrix of the coefficients is 


t Qi, %%2 - Q1°+* Sg,-1) 
Qi 1 2 - . . 


Qi1 Qie Qio 1 
0 (A-2) 


Si¢g,-1) 
Sie Gg-1) * Si(g,-1) 


fo = ae J;) and the right-hand sides of the equations contain terms 
J5-j k-1 


y Pi(j +k) tig +k) | J Qj +h)? corresponding to the expected duration of the look that discovers 


the object, which we shall ignore. The inverse of this matrix is easily shown to be 
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whence, corresponding to (5), we have 
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The first of these three expressions for Tij can be used in all cases by defining Qio = Qiy =0, 
i 


which can be interpreted to mean that the object, if it appeared before the beginning of the 
search period, would certainly be detected by or at the beginning of the search, and that, if it 
appears during the search, it is certainly detected by or at the end of the search. 

Substitution of (A-4) into (2) determines the Lagrange multipliers L;; whence T,; can 
be expressed explicitly as in (6). We can substitute the resulting Ty into (A-1), obtaining an 
expression analogous to (7), to be minimized by an appropriate choice of the J i subject to the 
condition analogous to (3) 
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However, the carrying out of this minimization presents difficulties (unless the tj are 
periodic as functions of j). 

One may further ask what set of tis satisfying (A-5) minimizes the expected delay. It 
would be desirable to know under what conditions the optimum set of ty is independent of j. 
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A GENERAL THEORY OF MEASUREMENT 
APPLICATIONS TO UTILITY* 


Johann Pfanzagl! 





For a long time measurement was confined to sciences such as geometry, astronomy, 
and physics. It is, therefore, quite natural that the theory of measurement was restricted to 
the special circumstances usually encountered in these sciences, to additive magnitudes such 
as "mass" or "electrical resistance" or "length." In psycho-physics, psychology, or welfare 
economics, on the other hand, additive magnitudes are hardly present. Therefore, new methods 
had to be developed—in psycho-physics, for example, the method of "bisection," successfully 
applied by S. S. Stevens and others (e.g., [16, 17]) in obtaining scales for subjective magnitudes 
such as pitch, loudness, etc., and in econometrics the method of Morgenstern-von Neumann 
for measuring subjective utility. 

Because these methods are not covered by the traditional theory of measurement, it 
was argued (especially in discussions concerning measurement in psycho-physics and psy- 
chology, e.g., [3]), that they do not lead to genuine measurement but only to something in the 
nature of an ordinal scale. In the field of psychology, especially, this criticism was supporied 
by the fact that there was no adequate theory as well founded as the theory of measuring addi- 
tive magnitudes. The situation concerning the measurement of utility is quite the opposite of 
course: Here, we have an incontestable theory, developed by Morgenstern and von Neumann 
[11], and it is rather the empirical investigations which are lacking. 

In the following sections we give a brief outline of a general theory of measurement 
and then its applications to the measurement of utility. 


OUTLINE OF A GENERAL THEORY OF MEASUREMENT 
The general aim of "measurement" is to map a set M on the set of real numbers in 

such a way, that—to the greatest possible extent—conclusions concerning the relations between 
elements of M can be drawn from corresponding relations between their assigned numbers. A 
trivial example: If a set is ordered, the mapping is performed in such a way that the order- 
relation between the elements of M is reflected by the order-relation of their assigned num- 
bers, or if there is an additive operation defined between each pair of elements of M, the 
mapping is performed in such a way, that the assigned number of a sum of elements in M 
equals the sum of the assigned numbers of these elements. 
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We now consider an ordered set M, for which there is defined an operation, called a 
"metric" operation, which fulfills the following set of axioms: 

1, Existence Axiom: For each pair of elements a,b¢€ M there exists a unique 
element ao be M. 





< 


< 
2. Monotonicity Axiom: If a rd a', then aob me a'ob for all be M. 





< < 
If b . b', then aob ~ aob' for all ae M. 
> 


3. Continuity Axiom:! The operation aob is continuous for both a and b. 





4. Bisymmetry Axiom: (aob) 0 (cod) ~ (aoc) o (bod). 





THEOREM 1. An ordered, connected set M, for which a metric operation is defined, 
can be mapped into the set R of real numbers (a € M, a ~ a* € R) in such a way that: 
1, The mapping a — a* is continuous, 


< 
2. The mapping a — a* is monotone, i.e., a S b implies a* 


vil A 


b*, 


3. The operation 'o" is mapped isomorphically on a linear operation: 
(aob)* = pa* + qb* +r. 


Mappings of this kind are unique up to linear transformations. An exact proof of 
Theorem 1 is given in [13, pp. 49-51]. If one does not insist on ultimate mathematical 
precision, one can assume in advance the possibility of a preliminary mapping a— a'e R 
which is continuous and monotone. Then, by (aob)' = F(a',b') a function is defined which is - 
due to the metric axioms - unique, continuous, monotone and bisymmetric: 


F[F(a',b'), F(c',d')] = F[F(a',c'), F(b', d')]. 


In [1] Aczél has proved that functions of this kind can be expressed by a function f(x") in the 
following way: 


F(a',b') = £7) [pf(a") + af(b') + r]. 
If we put x* = f(x') it follows immediately, that: 
(aob)* = f[F(a',b')] = pf(a') + qf(b') + r = pa* + qb* +r. 
As the scale is unique up to linear transformations, the magnitudes p and q are uniquely 
determined by the operation "o."" They are invariant under any linear transformation. 
(According to the monotonicity of the metric operation both p and q have to be positive.) In 


the general case p + q ~ 1, r is not invariant under linear transformations. Therefore, if the 
origin of the scale is chosen in an appropriate manner, r vanishes, so that (xoy)* =px* +qy*. 
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For a precise mathematical statement of this Axiom see [13, p. 20]. 
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Hence, in the regular case (p + q ~ 1), there exists a natural origin, and the scale is 
unique up to multiplication by a constant. 

In the singular case p + q = 1, even the constant r is invariant under linear transfor- 
mations and therefore uniquely determined by the metric operation "o."' Hence, no natural 
origin is fixed by the metric operation itself in the singular case: additional aspects have to be 
taken into account in order to fix the origin and to achieve a scale unique up to linear transfor- 
mations. (For instance, if M has a smallest element, it is quite natural to require that this 
smallest element shall be assigned the number 0.) 

A very important example of the singular case is the reflexive operation, for which 
aoa ~a for each ae M. The fact that p + q= 1 in this case follows immediately from relation 


a* = (aoa)* = pa* + qa* +r = (p + qga* +r. 


Furthermore, we conclude from this relation that r = 0 for all reflexive operations. 
The question as to under which circumstances two different metric operations lead to 
the same scale is answered by the following, Theorem 2: 


THEOREM 2: Two different metric operations "to" and "e" lead to scales identical up 
to linear transformations, if for each quadruple a, b, c, de M the following isometry-relation 
holds: (aob) @ (cod) ~ (aec) o (bed). 


For a proof of Theorem 2 see [13, p. 24]. The precise meaning of Theorem 2 is the 
following. Each of the operations "o" and "e" leads to a scale, which will be designated by 


a—a* and a—a** respectively: 


(aob)* 


P,2* + G,b* + Fe» 


(aeb)** pya** + qyb** + 1r,. 

Theorem 2 states that the two scales a* and a** differ only by a linear transformation if the 
two operations are isometric. The parameters Po» P43 Ag» 945 and To Ty will, in general, be 
different, of course. If at least one of the two operations is non-singular(e.g., Pota*)), we can 
choose a scale such that ro* 0: (aob)* = P,a* + q,*. If the two operations are isometric, 
then also r= 0: (aeb)** = p,a** + q,b**. This means that the natural origins are identical 
for both scales. 

The general theory outlined above covers, of course, the traditional case of an additive 
magnitude. From associativity and commutativity, the bisymmetry follows immediately, so that 
each additive operation is a metric operation in the sense defined above. Therefore, according 
to Theorem 1, a continuous and monotone mapping exists, such that (aob)* = pa* + qb* +r. 
Furthermore, as the additive operation is associative and commutative, p = q = 1 holds. 
Therefore, the origin of the scale can be chosen such that r = 0, so that we finally obtain a 
e } scale, for which (aob)* = a* + b* which corresponds to the traditional result. 

Also bisection, widely used by S. S. Stevens and others (e.g., [16, 17]) to construct 

} subjective scales of loudness, pitch, etc., seems to be a metric operation. The only relation 
whose validity could be questioned is the bisymmetry. The approach in the experiments was a 
purely empirical one, and there was no precise knowledge as to what conditions the operation 
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of bisection had to fulfill in order to allow for the construction of a cardinal scale. Therefore, 


no information concerning the validity of the bisymmetry axiom exists. But there is at least 
an indirect indication that the bisymmetry axiom might be valid. In the case of bisection, 


a hand ee an 





a aob 


aob is the magnitude midway between a and b, ao(aob) the magnitude midway between a and 
(aob), (aobjob the magnitude midway between (aob) and b. As the bisection is also com- 
mutative and reflexive, it follows, from bisymmetry, that (aob) should be midway between 
ao(aob) and (aob)job. According to the experiments performed, this relation seems to hold at 
least for the measurement of pitch. This is a hint that the bisymmetry axiom holds in the case 
of pitch and that bisection is therefore a metric operation.” In this case, according to our 
Theorem 1, a mapping exists, for which (aob)* = pq* + qb* + r. As stated above, the bisection 
is reflexive. It is therefore a singular metric operation, and no natural origin is determined 
by the operation itself. The origin can, nevertheless, be assigned in quite a natural way by 
assigning the number zero to the most extreme sensation at the lower bound (the least degree 
of loudness, the deepest pitch, etc.). However, the origin might be fixed arbitrarily; since the 
operation is singular, this has no influence over r: According to the assumption of reflexivity, 
r= 0 in any case. Together with commutativity, reflexivity yields that p = q = . Hence, there 


exists a scale, unique up to linear transformations, such that 


(aob)* = — (a* + b*). 


wl 


This is exactly the result Stevens assumed to hold. 
An additional example of the application of the general theory of measurement is the 
measurement of utility, which will be treated in detail in the next section. 


THE MEASUREMENT OF SUBJECTIVE UTILITY ACCORDING 
TO MORGENSTERN AND von NEUMANN 

In the following remarks it will be shown that the general theory outlined above can be 
applied to the measurement of subjective utility according to the concept of Morgenstern and 
von Neumann. 

M is interpreted as the set of the situations to be valued. With Morgenstern and von 
Neumann we assume that a complete order is defined for the elements of M; (see [11, p. 26, 
axiom 3A]). 

In [13, pp. 53, 54] it is shown that an ordered set M which fulfills the axioms of 
Morgenstern and von Neumann is connected and that the operation "agb" fulfills for each 
ae(0,1) the metric axioms. Therefore, the operation signified by a can be regarded as a 
special case of the metric operation"o." Furthermore, our isomeiry-relation (aab) p(cad) ~ 
~ (afc) a(bfd) is fulfilled for any a, 8 «(0, 1). This has the following consequences: 





2On the other hand, experiments reported by Gage [6] suggested tha‘! the relation may not hold 
for loudness. See also [12]. 
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1. For the construction of a cardinal scale it is sufficient to consider alternatives with 
any fixed probability a which can be chosen at will. It is not necessary to use alternatives 
with differing probabilities. For the special case of alternatives with the subjective proba- 
bility - this has already been pointed out in 1926 by Ramsey [14]. As the operation "ac b" is 


reflexive, we get a cardinal scale of utility U(x), for which U(aa b) = pU(a) + (1-p) U(b). This 
scale is unique up to linear transformations. 

2. As the isometry-relation is fulfilled we get, on the basis of alternatives with a 
probability 8 =a, identical scales (i.e., scales differing only by a linear transformation) of 
utility U(x), but different weights p', (1-p'). Therefore, U(agb) = p'U(a) + (1-p") U(b). 

Our assumptions suggest that the weights p and p" should be interpreted as subjective 
probabilities which are assigned to the objective probabilities a and 8; respectively, p = s(a), 
p' = s(8). The assumption that the operation aab fulfills not only the metric axioms but all 
axioms of Morgenstern and von Neumann, implies the identity of the subjective probability 
s(a@) with the objective probability a. Axiom 3:C:b [11, p. 26] states that (aab) Bb ~ aaBb. 
This implies s(a8) = s(a) - s(8). Together with s(1- a) = 1-s(a@), this leads to the result: 
s(a@) =a. 

However, it is not absolutely necessary to suppose that the axioms of Morgenstern and 
von Neumann are fulfilled. If our (weaker) metric axioms are fulfilled, this is sufficient for 
the construction of a cardinal scale of utility of the following form: U(aab) = s( a) U(a) 

+ (1-s(a@)) U(b). If, in addition, the isometry-relation is also fulfilled, then the same scale of 
utility is obtained, whatever probability is taken as a base for the alternative aa b. 

These assumptions are weaker than those of the system of axioms of Morgenstern and 
von Neumann. In spite of that, they permit us to derive all relevant results concerning the 
scale of utility, admitting, however, a divergence between subjective and objective probability. 

It should be noted that, on the basis of the general theory of measurement, the sub- 
jective probability assigned to a arises quite naturally out of the procedure of measurement. 
It is not necessary to anticipate its existence and other properties in a separate system of 
axioms in addition to the axioms concerning the operation aab. It would be worthwhile to 
undertake a critical examination of the consequences resulting from this fact for the theories 
which—following the example of Ramsey—start from a common axiom system for utility and 
subjective probability, such as those of Savage [15] and Luce [8]. 

Until now, the validity of the metric axioms and of the isometry-relation is derived 
from the axioms of Morgenstern and von Neumann. We will now try to check the evidence of 
the metric axioms and of the isometry-relation directly. We shall restrict this undertaking to 
measuring subjective utility with regard to different quantities of an identical commodity, e.g., 
of money. By this restriction the problem can be treated more precisely and the analysis of 
the necessary assumptions can be carried through more thoroughly. 

In the subsequent analysis the elements x eM therefore stand for different quantities 
of the same commodity. Furthermore, an event is given, the occurrence of which is uncertain. 
Let us designate the occurrence of this event by P, its non-occurrence by P. The wager 
"aPb" with a,be M means: The individual in question gets quantity a, if P occurs, quantity b, 
if P occurs. Let {P,P} be an event which can be repeated indefinitely, e.g., a random experi- 
ment. (In principle the repetition could also consist in the fact that a unique event {P,P} is 
judged by different persons with identical preference scales.) Let m(aPb) be the quantity 
which has the same subjective utility as the wager "aPb". As is shown even by every-day 
experience, the quantity m(aPb) is not uniquely determined without further conventions. 
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In order to determine this quantity uniquely the following method is used (see [10]). Let 
w(x,y) be the probability that the individual prefers x, if faced with the alternatives x and y. 
m(aPb) is then defined by w(m(aPb),aPb) = a. Of course, any deviation of m(aPb) from the 


objective expected value (in the sense of the probability calculus) of the wager aPb is possible. 

After these preliminary remarks we shall now proceed to an examination of the evidence 
of our metric axioms. 

1. Existence axiom: To the two elements a,b € M we assign the element m(aPb) « M. 
This merely assumes that the above outlined definition of m(aPb) is meaningful. 

2. Monotonicity axiom: The validity of the monotonicity axiom seems to be evident 
since we have confined the definition of m(aPb) to a,b being different quantities of the same 
commodity. 

3. Continuity axiom: Also this axiom seems evident in consequence of the restriction 
to different quantities of the same commodity. 

4. Bisymmetry axiom: m[m(aPb) Pm(cPd)]= m[m(aPc) Pm(bPd)]. Let P and Q be 
two different events. Then certainly (aPb) Q(cPd) = (aQc) P(bQd). For, whichever of the 
combinations PQ, PQ, PQ, PQ is realized, the result will always be the same for both 
(aPb) Q(cPd) and (aQc) P(bQd), namely: PQ—a, PQ—b, PQ —c, PQ—d. In view of this fact, we 
have used the symbol of identity "=" in order to distinguish this relation from the equivalence 
"~" defined in terms of utility. The decisive assumption is that uPv ~m(uPv) and xPy ~ m(xPy) 
implies (uPv) Q(xPy) ~m(uPv) Qm(xPy). By means of this relation we can deduce from the 
identity stated above that m(aPb) Qm(cPd) ~ m(aQc) Pm(bQd). Furthermore, we must assume 
that uPv ~xQy (i.e., w(uPv, xQy) = —) implies m(uPv) = m(xQy). Then we get m[m(aPb) Q 


m(cPd)] = m[ m(aQc) Pm(bQd)]. This isthe isometry-relation which guarantees that the events 
{P,P} and {Q, Q} lead to the same scale of utility. To obtain the bisymmetry axiom we need 
only assume that {Q, Q} is a repetition of the experiment {P,P}, independent of the result 
which this experiment has yielded. Then m(xQy) = m(xPy), and thus m[m(aPb) Pm(cPd)] 
= m{m(aPc) Pm(bPd)]. 

5. Reflexivity axiom: m(aPa) = a seems evident, as the wager aPa leads with 




















certainty to the result a. 

The relation m(aPb) = m(bPa) does not involve an additional assumption; it holds a 
priori as a consequence of the designation we are using here. 

Summing up, we can say that the validity of the metric-axioms, plausible as they are, 
still needs a critical empirical verification. The possibility does not seem to be excluded that 
e.g., experiments with wagers xP*y with a subjective probability 5 will show that m(xP*y) = 
= p Max(x,y) + (1-p) Min(x,y); p being any number between 0 and 1. For p # B the bisymmetry 
axiom is not fulfilled; in this case, the construction of a metric-scale of utility on the basis of 
the wagers xP*y would therefore be impossible. In spite of that, the behavior of the person in 
question could not be called irrational. 


THE CONSISTENCY AXIOM 
The following section considers the meaning and the implications of an additional 
assumption concerning the evaluation of wagers 


Consistency axiom: m[(a+ c) P(b + c)] = m(aPb) + c. 
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The axioms dealt with above were exclusively concerned with the evaluation of alter- 
native events. In the case of the consistency axiom the evaluation of conjunctive events is 
considered. We designate the conjunctive connection of two events x,y by "x&y" (both x and 

y). Then certainly (a + c) P(b + c)=(aPb) & c. The right-hand side of the equation means: 

the amount c is paid in any case, and, in addition, the amount a is paid in the case of P, and 

the amount b in the case of P. This is identical with paying the amount (a +c) in the case of 

P, the amount (b +c) in the case of P. If one assumes that aPb ~ m(aPb) implies 

(aPb) & c ~ m (aPb) & c, then one gets: (a + c) P(b +c) ~ m(aPc) & c = m(aPc) +c. This 

means by definition: m[(a + c) P(b + c)] = m(aPc) +c, as stated by the consistency axiom. 

Though the consistency axiom was up to now not stated explicitly anywhere, it can be 
shown that several authors are really tacitly assuming it. In the following, we will show this 
has been so in the case of Friedman and Savage [5], Mosteller and Nogee [10], and von Neumann 
and Morgenstern [11]. 

Using the designation adopted in this paper,>? a statement by Friedman and Savage in 
[5, p. 290] reads: "If m(a'ab') is greater than x, the consumer unit (purchaser) prefers this 
particular risk (namely, the participation in the lottery a'ab') to a certain income of the same 
actuarial value and would be willing to pay a maximum of [m(a'ab') - x] for the privilege of 
'gambling'."' This statement by Friedman and Savage is confined to the case where x is the 
actuarial value of the lottery (a'ab'). But this argument, if valid for the actuarial value, is 
obviously valid for any value x, so that we get the equivalent statement: "If m(a'ab') is 
greater than x, the purchaser prefers this particular risk (namely, the participation in the 
lottery a'ab") to a certain income of the value x and would be willing to pay a maximum of 
[m(a'ab') - x] for the privilege of 'gambling'." 

If the purchaser pays an amount [m(a'ab') - x] for the lottery ticket (a'ab'), he 
actually faces the chances of getting the amount [a' - m(a'ab') + x] with probability a and 
[b' - m(a'ab') + x] with probability (1-a@). The utility of this risk is—according to the quota- 
tion above—equal to the utility of the amount x: The purchaser is willing to pay [m(a'ab') - x] 
for the lottery ticket; therefore: U ([a' - m(a'ab') + x] a [b' - m(a'ab') + x]) 2 U(x). But the 
purchaser is not willing to pay more than [m(a'ab') - x]. Therefore: U ([a' - m(a'ab') + x] 

a [b' - m(a'ab') + x]) = U(x). It follows: U ([a' - m(a'ab') + x] a [b' - m(a'ab') + x]) = 
U(x). This means: x = m (|a' - m(a'ab') + x] a [b' - m(a'ab') + x]). By putting m(a'ab')-x = 
= c we obtain: m(a'ab') - c = m[(a'-c) @ (b'-c)]. If we put further a'=a+c, b'=b+c, we 
obtain: m[(a +c) a (b+ c)]= m(aab) +c. Therefore the validity of the argument used by 
Friedman and Savage implies the validity of the consistency axiom. (As a and b are incomes 

in the case discussed by Friedman and Savage, we have assumed in the derivation given above, 
that there is no other amount of money to be taken into account. But this was only for sake of 
brevity. The same argument holds if one regards s + a'-m(a'ab') + x, s + b'm(a'ab') + x and 

Ss +x instead of a'-m(a'ab') + x, b'-m(a'ab') + x and x respectively. 

Mosteller and Nogee [10, p. 399], in discussing the results and conditions of their 
experiment, consider also the "effect of the amount of money in front of the subject upon his 
decisions" to take part in the gamble or not. One of their statements is this: "One possible 
criticism could be that the subject changes his utility curve with these changes in capital, so 
that each decision he makes depends on the amount of money he has on hand at that particular 
moment." It seems remarkable that Mosteller and Nogee are thinking only of the possibility 





3The translation in the designation is: I, ~a', I, ~b', I*—~m(a'ab'), Ix. 
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that the available amount of money could change the utility function of the player. They do not 
consider the possibility that the willingness to accept a particular game could depend on the 
amount of money held by the individual, even if his utility function were unchanged. And they 
are only thinking of "the amount of money in front of the subject""—not of the amount in his 
pocket. This becomes especially clear when they are talking (in a different connection, p. 403) 
of one of the participants, who became unemployed during the course of the study. It appears 
that these authors had neglected the possibility that this might have changed his behavior con- 
cerning the participation in these games. If this were true, they would have assumed that the 
utility of a special game is independent of the money held by the individual (perhaps with an 
exception regarding the money immediately involved in playing). This assumption, however, 
is equivalent to the consistency axiom. 

The decisive question seems to be this: $10 a $1 is a lottery, offering an amount of 
$10 with probability @ and an amount of $1 with probability 1-a@. If an individual is willing to 
pay $2 for a ticket of this lottery (and not more than $2), this does not necessarily mean that, 
for this individual, $10 a $1 is equivalent to $2, i.e., m(10 a 1) = 2. It actually means that the 
status quo has the same utility as a lottery, which leads to the status "quo + $8" with proba- 
bility a and status "quo - $1" with probability 1-a@. (We must insert $8 and - $1 as net out- 
comes of the lottery, since we have to deduct the $2 spent for the lottery ticket from the gross 
prices of $10 and $1, respectively.) If we assume that for this purpose the status quo can be 
described essentially by stating the amount of money s held by the individual, then the willing- 
ness to pay at most $2 for the lottery ticket really means that m[(s + 8) a (s-1)]=s, not 
m(10 @ 1) = 2, as stated above. Both statements are equivalent only if the consistency axiom 
holds. 

It would require a thorough investigation as to what the real meaning of experiments 
like those performed by Mosteller and Nogee [10], or by the Applied Mathematics and Statis- 
tics Laboratory of Stanford University (e.g., [4] or [18]), could be, if the consistency axiom 
were not fulfilled. 

Finally, the consistency axiom is also used in the theory of games. A good example 
is von Neumann-Morgenstern's introduction of the concept of "strategic equivalence" in [11, 
p. 245 ff]. My attention was drawn by Oskar Morgenstern and Harlan D. Mills to the results 
obtained by Kemeny, DeLeeuw, Snell, and Thompson [7]. According to a quotation in [9, p. 72}, 
in examining the restrictions imposed on the utility function by the concept of strategic 
equivalence these authors have shown in [7] that the utility function-assuming strict mono- 
tonicity and differentiability must be one of the types stated in Theorem 3. 


THEOREM 3. If utility is measurable by a cardinal scale such that U (m(aPb) ) = pU(a) 
+ (1-p) U(b), if furthermore the consistency axiom holds, and if the utility function is assumed 
to be continuous and strictly monotone, then the function is of one of the following types: 


U,(x) = Ax +B A>O, 


U,(x) = AA*~ +B A>0,A>1 or A<0,0<A<1. 
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If U(x) is standardized such that U(0) = 0,U(1) = 1, we obtain: 





U,(*) =x 
x 
U,(x) = tA ‘; 
1-A 


In U, (x) we must distinguish between the case A > 1 and 0<A< 1. In both cases U,(x) is 
monotone increasing; however, the marginal utility U'(x) decreases only in the case O<A<1 
with increasing x. In this case 





lim U,(x) = 1 
xz? @ 1-A 


These statements concerning the shape of the scale of utility cannot be verified imme- 
diately. They can, however, be put in a way which allows an experimental proof or disproof. 
The whole theory of the measurement of utility depends on the fact that the magnitude m(aPb) 
can be clearly determined by means of experiments. This, however, permits us to transform 
the statements concerning the shape of the scale of utility into statements concerning directly 
observable magnitudes: It follows from the two solutions U,(*), U(x), that m(aPb) must be a 
function of one of the following types: 


pa + (1-p)b 


ul} 


m,,(aPb) 


m,(aPb) = log, [pa® + (1-p) Ned . 


Both functions fulfill, as can be easily verified, the metric-axioms and the consistency axiom. 
Furthermore the isometry-relation is fulfilled if 4 is independent of the special event {P,P} . 
Given the subjective probability p, the behavior of the individual can be described by only one 
parameter, A. 
Proof of theorem 3: From U(m(xPy)) =pU(x) +(1-p)U(y) and m[(x+z)P(y +z) |= m(xPy) +z 
together we obtain the functional equation: 
pU(x + z) + (1-p)U(y + z) = u{ u7 {pus +(1-p)U(y)]+z}, O<p<1. 

We are considering continuous and strictly monotone solutions of this functional equation, 

Let y = 0: pU(x + z) + (1-p)U(z) = u{u7 pus) + (1-p)U(0)] + z} 

Let x = 0: pU(z) + (1-p)U(y + z) = u{u7 fpu(o) + (1-p)U(y)] +z}. 


By adding these two equations and using the functional equation stated above again, we obtain: 


u{u~ [pu(x) + (1-p)U(y)] + z} + U(z) = 


(1) 


= u{u~|[pu(x) + (1-p)U(0)] + z} + U{U~ [pU(O) + (1-p)U(y)] + z} 
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If we introduce the following: 
p[U(x) - U(0)]= é sis 
(1-p) [U(y) - U(0)] =n inv 
u{u-*[¢ + W(0)] + z} - U(z) = {(€,z), 
then it follows from (1) that 
£(E +n,z)=£(& ,z) +f£(n,z). 
{(€ ,z) is for each z a continuous and strictly monotone function of ¢ . 
Hence: £(€ ,2z)=€ (2), Lew uf u7 [e+ U(0)] + z} - U(z) = ¢ gz). We designate: | f 
U(x) - U(0) = W(x) and obtain: W{W"* (¢ ) +z} = ¢ v(z)+¥(z). For € = ¥(x) it follows: 
(2) U(x + z) = (x) g(z) + ¥(z). 


As ¥(x +z) = ¥(z + x), we obtain: (3 








W(x) p(z) + ¥(z) = ¥(z) p(x) + (x). 
[4 
As U(x) is strictly monotone, ¥(x) = 0 for x ~ 0, hence: 
P(z)-1_ 9-1 5, PM-1_,, q [6 
¥(z) ¥ (x) ¥ (x) 
This result introduced in (2) gives: [6 
W(x + z) = c W(x) W(z) + W(x) + W(z). } [7 
This is a functional equation of the type W(x + z) = F[W(x), ¥(z)]. The continuous and strictly 
monotone solutions of such a functional equation are essentially unique: If ¥ o(*) is a con- ( 
tinuous and strictly monotone solution, then all functions ¥ (ax) are solutions and only these 
(cf. Aczél (2, p. 120, 121]). We have to distinguish between two cases: 4 (9 
[10 
c= 0: general solution: ¥ o(*) = ax c 
, [11 
c = 0: general solution: ¥ (x) a= (Q*-1). 
c 
[12 
From W(x) = U(x) - U(0) we obtain the theorem to be proved. 
CONCLUSIONS ‘ [13 


The applications of the results of a general theory of measurement, developed in [13], 
to problems of utility show that the valuation of wagers with a constant probability is sufficient j 












for the construction of a cardinal scale. This cardinal scale allows for a divergence between 
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subjective and objective probability. Even this "weak" cardinal scale, together with the con- 
sistency axiom introduced in this paper, restricts the utility functions essentially to a one- 
parameter family. If, on the other hand, the consistency axiom is not valid, a more thorough 
investigation would be needed in order to determine the real meaning of the utility functions 
obtained by experiments with gambling. 


(1] 
[2] 


[3] 


[4] 
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[7] 
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ON THE SOLUTION OF PROBABILISTIC LEONTIEF SYSTEMS* 


Richard E. Quanat! 


Princeton University 


INTRODUCTION 

The objective of input-output analysis is to describe and analyze the interdependence 
among the sectors of the economy. The fundamental tool of input-output analysis is a system 
of simultaneous linear equations (called the "Leontief system" after Wassily W. Leontief) 
given by 


(1) (I- A)X= Y 


where I is the identity matrix, A the matrix of input coefficients, X the column vector of 
gross outputs in constant dollar terms in each of n industries, and Y the column vector of 
final demands. The input coefficient aij is the fraction of a dollar's worth of commodity i 
necessary to sustain the production of one dollar's worth of commodity j. Equation (1) thus 
states that the aggregate amount of commodity i produced minus the amounts of commodity i 
used up in producing all commodities equals some residual Yj which can be assumed for 
simplicity's sake to accrue to consumers. 

If the coefficients of A are known, and if one stipulates a final demand vector Y*, one 
can determine the gross outputs, which are necessary in order to permit consumers to con- 
sume the quantities indicated by the elements of Y*, by solving (1): 


(2) x = (1- A)7ly*. 


In general it is reasonable to assume that an observed aij is a random drawing from a 
probability distribution defined over aij. In this case the stochastic properties of the solution 
vector X become of interest since the probability distributions defined on the ai; implicitly 
define a probability distribution on X, provided that I - A is non-singular. It is relatively 
simple to construct confidence intervals for the solution if it is assumed that errors in the 
coefficients are normally distributed [1]. Given an estimate of the variances of the errors in 
the coefficients a,,, one can determine the variances of the elements of the solution vector [3] 
and establish weak confidence limits even if the distribution of the errors is not known [4]. 





*Manuscript received November 26, 1958. 
The author acknowledges his indebtedness to Professor Frank Anscombe for comments and to 


Dr. Fred Brooks for programming the sampling experiment described in Section 2 on the Mark 
IV computer. 
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However, in order to get relatively narrow confidence intervals, it is obviously not sufficient 
to rely on procedures such as Chebyshev's Inequality. It then becomes necessary to know the 
distribution of the solution. 

In certain special cases it may be possible analytically to derive the distribution of the 
solution, but in general this does not appear to be possible [5,6]. For this reason it was 
attempted to investigate the properties of the solution of a Leontief system by resorting toa 
large number of sampling experiments. The purpose of the present paper is to describe these 
sampling experiments and to discuss their results and implications. 


THE PROCEDURE FOR THE EXPERIMENTS 

In general terms the following procedure was observed. An initial hypothetical 3 x 3 
input-output matrix was specified. This matrix was considered to be the true matrix throughout 
the experiment. A bill of goods was also specified and was held constant throughout the exper- 
iments.! The magnitude and probability of errors in the coefficients were subsequently decided 
upon; i.e., discrete probability distributions were specified for the coefficients. The coefficients 
were then altered (i.e., perturbed or shocked) in accordance with the specified probabilities. 
A large number of random Leontief matrices was obtained in this fashion. Each one-of the 
matrices was inverted and for each one of these random systems of equations the solution was 
obtained. Finally, the relevant moments of the distributions of the inverse elements and of the 
solution were obtained. 


The Initial Matrix 
The initial true matrix was assumed to be 


[ .50 


15 
10 


and hence 


00 8 -.30 
-.15 90 
|--10 =-.05 








The bill of goods (assumed constant throughout the experiment) was assumed to be the vector 
(.03, .07, .10). 


The Nature of Errors 

The procedure of shocking A (or I - A) has been carried out for eleven different 
distributions of errors. These distributions are described in Table 1. It was assumed 
throughout that errors in coefficients are independently distributed. This assumption prob- 
ably causes some loss of realism in the experiments since it is often conjectured that errors 
in Leontief systems are mostly of the compensating type. The assumption was made none- 
theless in order to avoid certain programming difficulties for the Mark IV computer. 





I The bill of goods has no intrinsic significance and was chosen merely for the purpose of allow- 
ing the system to be solved. It is arbitrary in every respect. 
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TABLE 1 
Distribution Errors Assigned Probability of 
toI-A Errors 

I 05 O -.05 32 34 34 

I 05 O -.05 -59 -30 | 
Il 05 O -.05 -74 -20 .06 
IV 05 O -.05 -09 -30 -61 
V 05 O -.05 -04 -20 -76 
VI 05 O -.05 44 10 -46 
vil 05 O -.05 19 -60 21 
vill 02 O -.05 32 .34 .34 
Ix 05 O -.02 232 34 -34 
x .05 .02 0 32 34 34 
XI 0 02 0.05 34 32 34 

















The distribution of the errors ei were specified in such a fashion that aij + ei 20 for 


all i and j and = (aj, + ej) <1 for all j. This guarantees that the determinant |I - A - E|is 
i 


shielded away from zero, and no nonsingular system may be generated in the process of shock- 
ing matrices [7]. 


Sample Size 

For each distribution a sample of 100 matrices was obtained. Since each matrix has 
nine elements and each element can have one of three values, for any one of the error dis- 
tributions the population of matrices consists of 39 = 19,683 matrices. Each sample com- 
prises about 0.5 percent of the population and the resulting variances can be considered to be 
estimates of the population variance without the customary adjustment for finite population. 


The Method of Perturbation 

The probabilities assigned to the errors in the various distributions specify two 2-digit 
numbers a and b which have the property that randomly chosen 2-digit numbers x fall into 
the ranges 00 £=x<a, a=x<b, b =x =99 with the same probabilities as are assigned to the 
particular errors. The computer was instructed 

(1) to produce a two-digit random number, 

(2) to inspect the random number and augment the first element of I - A by one of the 
three possible errors depending on whether the random number falls into the first, 
second, or third range, 

(3) to move to the second element of I - A (proceeding rowwise), and repeat steps (1) 
and (2) until all elements of I - A have been shocked, 

(4) to repeat steps (1), (2), and (3) on the initial matrix until a total of 100 matrices 
have been generated, 

(5) to invert each matrix and to obtain the solution of each system of equations, 

(6) to calculate some relevant moments of the distribution of the solution. 
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Instead of random numbers the computer utilized pseudo-random numbers generated 
according to a scheme suggested by Votaw and Rafferty [8]. This procedure has the advantage 
that random numbers can be produced easily on the computer and that the resulting sequence 
of numbers exhibits randomness in the statistical sense for relatively short sequences. The 
length of the sequence utilized here was substantially less than the maximum allowable length. 


INTERPRETATION OF THE RESULTS 


TABLE 2 








Distribution | 0” a, a, Distribution | 07 





Sy 


a4 








eee 


.001653 
-001174 
-000844 
.001074 
-000704 
-002253 
-001003 
-000866 
-0008 52 
-000416 
-000430 


-11090 
-.95341 
-1.64732 
1.03270 
1.74996 
.08653 
-17327 
-.44033 
-52266 
-28313 
-.22769 


1.51694 
2.66106 
4.63685 
2.89739 
5.15147 
1.11278 
2.50098 
1.47235 
1.56553 
1.54318 
1.47215 





Mu Rods<2hne 


-001689 
.000526 
-000274 
.002794 
-003201 
-002477 
-000860 
-001506 
.000569 
-000159 
.001111 


-5852 
-8019 
1.2290 
-1651 
-.1636 
-6324 
-4542 
-7681 
2481 
.0547 
-5406 


3.5508 
3.5925 
5.1428 
2.6987 
2.1943 
2.7765 
3.5540 
2.9955 
2.2187 
4.0000 
2.6737 






































Table 2 contains the variances and the third and fourth standard moments of the error 
distributions; Table 3 contains the same moments for the distributions of the first solution 
element x,. The third standard moment, denoted by a3, is a common measure of skewness, 
although not necessarily an always relaible one since certain asymmetric distributions may 
have zero third moment around the mean [9]. The fourth standard moment, denoted a,, is a 
common measure of kurtosis, although again not necessarily a relaible one [10]. 

In order to examine the relationship between the moments in Table 2 and the moments 
of the distributions of the solution, the relevant moments of the eleven distributions were 
ranked, and rank correlation coefficients were obtained between the ranks of moments of the 
error distributions and the ranks of the moments of the distributions of the solution. 
Spearman's p has been employed throughout for this purpose [11].2 The notation employed 
can be illustrated by p(a 3e? “3x,) which denotes the rank correlation coefficient between the 


ranking of the error distributions according to their third standard moments and the ranking 
of the distributions of Xy according to their third standard moments. 





2Other coefficients of rank correlation are also available but Spearman's phas been preferred 
on account of its flexibility and ease of computation. p is defined for the interval from -1 to 1. 
Perfect agreement between the rankings is indicated by P= 1, perfect disagreement by Pp = -1. 
It would have been possible to examine the relationship between the moments of the original 
error distribution and the moments of the distributions of the solution by calculating the 
ordinary product-moment correlation. This procedure was discarded on the grounds that one 
would have to make unjustifiably strong assumptions about the relationship, i.e., assumptions 
of linearity or some particular type of curvilinearity. 
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The first test was intended to examine the possibility of some correlation between 
o » Ae, and a 4e* If such correlation existed, one could not be certain whether a high rank 


correlation between, say, @ 3e and As. might not be the result of a high correlation between 
1 


Ago and 0 on the one hand and between o2 


“ ° and e, on the other. The relevant rank cor- 


relations are 


2 " 

p(s, ag.) = -.1636 
p(o*, a,.) = -.2000 
e? 4e ‘ 
p(at,, @4,) = .2727 


In no case is it possible to reject the null hypothesis that Pp = 0.3 
Next it is desired to correlate the moments 0 1%, Ag, with the moments of the 


distributions of X41» Xgy and Xq- If, employing the same criterion of ranking, the distributions 
of Xo and Xq are ranked essentially identically with the distributions of X1) it is sufficient to 
correlate the moments 02 » Fo, By, with the moments of xy only. The degree of agreement 


among three or more rankings is measured by the coefficient of concordance W defined by 


12S 


Ws... 
m? (n° - n) 


where m is the number of judges (3 in the present case), n the number of rankings (see foot- 
note 3), and S = Z(Z qj; - mM)? where qi; is the rank assigned to the ith item by the jth judge 
i j 


and M is the grand mean of the ranks. When n > 7, the quantity m(n - 1)W is distributed as 
x2 with n - 1 degrees of freedom [11, pp. 80-81]. The observed values or W are 


W(o2) = .9058, 
W(a.,) = .9324, 
W(a4,) = 9091. 
All three W's are significantly different from zero at the .001 level of significance. 


3Under the null hypothesis the quantity 


4 


p e-8 ‘ 
(1 - p>) 


is distributed as t with n - 2 degrees of freedom [11, pp. 47-48]. 
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Ranking the eleven error distributions and 
the resulting distributions of Xy according to the 
size of their various moments and finding all pos- 
sible rank correlations one obtains the results in 
Table 4 the entries of which are the rank correla- 
tions between the error distributions and the 
distributions of Xy ranked according to the criteria 
shown in the first row and first column respectively. 

Only three of the rank correlation coef- 
ficients in Table 4 are significant at all. The 
coefficient in the first row and second column is 
significant on the .20 level, the one in the second row and second column on the .05 level and 
the one in the third row and second column on the .10 level. Only Pass . 4.) can really be 





























considered to be a significant result. It means that the skewness of the original distribution 
tends to be transmitted to the distribution of the solution.* It is somewhat surprising that the 
variance of the original distribution does not seem to be transmitted as faithfully to the dis- 
tribution of the solution. The obvious reason for this is that the infinite series I + A + a2 
—— rather slowly on account of the relatively large size of the elements aij. Since 
(I - A)” ns oo Oe tat ..- when the column sums of A all add up to less than unity, very 
rapid convergence would mean that the inverse of the Leontief matrix would be dominated by 
the first few terms. In such an instance one could reasonably expect that the larger the 
variances of the elements of A, the larger would be the variances of the inverse elements 
and hence of the solution.> The fairly large values of the P's in the second column and the 
first and third rows cannot be imputed substantive significance since the distributions of Xy 
ranked according to their various moments are fairly highly correlated with each other. The 
relevant rank correlations are lox, . 3x) = -.4318, Hox, 4 4x.) = -.7364, Plas, = 4x,) 
= 9182. 

Lognormal distributions were fitted to each of the eleven frequency distributions by the 
method of moments. The fitted distributions were compared with the observed ones on the 
basis of the criterion 


(e; - £)” 
ce 


where ej is the expected frequency in the ith class interval and fi the observed frequency, 
and where the quantity K is distributed approximately as 2,6 





4The rank correlation coefficient is negative because the error distributions are attached to 
I-A and not to A, 

5It is pee, 02 , however, that if one omits Distribution V and reranks the remaining distribu- 
tions, pio2 ? becomes .6485, which is significant at the .05 level. 


6No class pile. was allowed to stand by itself if the expected frequency in it was less than 4. 
Usually 5 is used as the critical number, but it is sometimes argued that this is a much too 


conservative criterion. Cochran [12] argues that one may allow to stand class intervals with 
expected frequencies as low as l. 
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TABLE 5 
Distribution Degrees of K Probability of a 
Freedom Worse Fit 

I 2 0.4631 0.794 

sf 1 1.8036 0.178 
I 3 3.4335 0.330 
IV 2 1.7185 0.424 
Vv 4 3.3922 0.495 
VI 2 1.6862 0.430 
vil 2 1.3831 0.501 
Vil i 1.7190 0.190 
x 1 2.6274 0.105 
x 1 1.2753 0.529 
XI a 2.7296 0.604 




















Table 5 gives for each distribution the number of degrees of freedom, the observed 
value of K, and the probability of obtaining a worse fit on the assumption that the null hypoth- 
esis is true. Averaging the probabilities of obtaining worse fits, one gets 0.416. This is quite 
close to the 0.500 that one would expect to observe in an indefinitely large number of such 
experiments if the null hypothesis were true. One may therefore conclude that the lognormal 
distribution is a good empirical description of the distribution of Xy- 

The discovery that the lognormal distribution can be utilized as a description of the 
solution is important for several reasons. First, it justifies the use of a, because in the case 
of the lognormal distribution the common alternative measure of skewness, namely (mean - 
mode)/standard deviation, is not a satisfactory one [13]. Secondly, one can explain why the 
rank correlations between "3x, and 0, %,,%4, on the one hand and between @ 4x, and 


0 7%, %4, On the other are roughly of the same order of magnitude: it can be shown that 


there is a functional relationship between a. and @ 4 for the lognormal distribution [13]. The 
third and most important consequence of this finding is that the critical values of the lognormal 
distribution can be utilized in order to establish approximate confidence limits for the solution. 
The lognormal distribution is a three-parameter family of distributions given by 


2 
2 x-a 
f(x) = Waa exp[ (1/2¢) (lee : ) ‘ 


Confidence limits can be established by 
x+tj0< & <x+tyo 


where x and Oo are the estimated solution and its standard deviation and where ty and ty are 
given by 








2 
- oc-& 
oc-%5 








-1 
t. = 
1 9 1/2 
c 
(.< - 1) 
2 
%c-F 
Js - 1 
9 = 
2 1/2 
c 
(ec. 1) 


where @¢ is the point on the normal distribution corresponding to the desired level of signifi- 
cance [14]. Of course, in order to utilize tables of the lognormal distribution, one must know 
the skewness of the distribution. An approximation to the skewness can be obtained by follow- 
ing the procedure of Etherington [3]. Consider the system of n simultaneous equations 


(3) BX = Y 


where B is non-singular. Let B be the true matrix and B + E the observed matrix where E 
is the matrix of errors. Let D, be the determinant of the matrix C defined as the matrix B, 
with the first column replaced by the vector Y (assumed to be free of error), let Do be the 
determinant of B, and Ej and Eo the corresponding errors in the observed determinants. 
Clearly, 


n n-1 n 2 Zt n 
(4) IB+El=(Bl+) Bi+)' ) Bis...4) D1... D) BoM oie, 
i=1 i=1 j=2 i=l j=2 k=n-1 
i<j i<j<...<k 


where pi: - ok is the determinant of the matrix obtained from B by replacing its ith column 
by the ith column of E, its jth column by the jth column of E, etc. Then 


n n-1 n 
(5) E,=) B+)’ )' sis... +e 
i=l i=l j=2 
i<j 
and 
n n-l n i. 
(6) E,=) ce Ps. 
i=2 i=2 j=3 
where C--- (i,j, ...,k = 2) is the determinant of C with its ith, jth,...,kth columns 


replaced by the corresponding columns of E. Expanding the determinants in (5) and (6) in 
terms of the columns containing the errors one has 
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where Bi is the cofactor in Bi corresponding to the (ri)th element, Bi ,8j 
Cc ij 
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the (sj)th element in B ij and similarly for ci 
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The error in the solution can be expressed approximately as 


ri?’ 


cil 


ri,sj ~ Cri sj’ otc. 


contains any errors and thus Bi =B 
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Substituting in (9) one obtains as a crude approximation 


n n 


n n 
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is the cofactor of 


etc. Clearly none of these cofactors 


Assuming that all ej are distributed independently of each other with variance of the various 


moments of A(x;) can be obtained by performing the appropriate integrations. One thus 
obtains (to a very rough approximation) 
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where ra denotes the third moment around the mean. Obviously (3-11) is the only crudest 
approximation. Fortunately such an approximation can be calculated to any order, although 
the complexity of the necessary algebra increases rapidly. 


CONCLUSION 

The two most important conclusions from a large number of sampling experiments with 
probabilistically generated Leontief matrices are (1) that the skewness of the errors in the 
Leontief matrices tends to be transmitted to the solution and (2) that the lognormal distribution 
provides a fairly adequate description of the distribution of the solution, irrespective of the 
distribution of the original errors. Since the variance and the third standard moment of the 
solution can be calculated to an approximation from the moments of the original error distri- 
bution, one can utilize the critical values of the lognormal distribution to establish confidence 
intervals for the solution. Such an interval estimate is probably preferable to a point estimate, 
even though its accuracy is impaired by (1) the fact that the lognormal distribution is merely 
an approximation to the exact distribution of the solution, (2) the fact that various approxima- 
tions have to be used in calculating the relevant moments of the solution, and (3) the fact that 
the use of the critical values of the lognormal distribution may be illegitimate unless the 
population values of o2 and a 3 are known. 
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RELATIONS AND GRAPHS 

In a recent series of papers by one of us [4, 5, 6], the concept of domination of 
imputations in n-person games as developed in the now classical book by von Neumann and 
Morgenstern [3] was studied. In this formulation, each imputation is replaced by an element 
a, b, c,... of some universal set P. If a dominates b, then we write a >b. In this setting, 
we do not have a>a for any element a in P. Thus the relation > is irreflexive. The subset 
S of P is a solution of the irreflexive relation if the following two conditions hold: 

(1) For any two elements a, b in S, it is not true that a > b. 

(2) For any ae P-S, there exists beS such that b> a. 
These conditions may be rephrased: 

(1) S is internally stable. 

(2) S is externally stable. 
An irreflexive relation is solvable if it has at least one solution. 

It is natural to propose the following two problems: 








Problem 1. Characterize solvable and uniquely solvable irreflexive relations. 
Problem 2, Find all the solutions of a solvable irreflexive relation. 


Our main purposes are to supply an algorithm for Problem 2 for the case of finite 
irreflexive relations and to develop some properties of a generalization of the concept of 
solution. The procedure for finding all the solutions is based on an observation concerning 
cliques [1]. 

In different terminology, the same subject matter is mentioned in the book on graph 
theory by Kénig [2]. A directed graph or digraph D consists of a collection P of points 
together with a prescribed subcollection of the set of all possible ordered pairs of distinct 
points. If (a, b) is in the digraph D, then this ordered couple is called a (directed) line of D 
and is written ab; we then say a is adjacent to b and b is adjacent from a. Kénig defines 
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a point basis of the second kind as a set of points of a digraph such that there is no line joining 
any two of these points, and any point not in this set is adjacent from a point in this set. He 
then poses the problem of finding for a given digraph a minimal point basis of the second kind 
(solution). 

From now on, we use the language of graphs since it is more concise for the present 
purpose than that of binary relations. Thus we refer to the solutions of a digraph, the solva- 
bility of a digraph, etc. In [3, 4, 5, 6] there is a succession of sufficient conditions for a 
digraph to be solvable. We recapitulate some of these, and add some additional observations. 
We conclude by showing that a solution of a digraph and a point basis (of the first kind) of a 
digraph are special cases of an r-basis. If n is the number of points of D, then a point basis 
is an (n-1)-basis, while a solution is a 1-basis. An indication that the concept of an r-basis 
may prove to be fruitful arises from the fact that there already exists in the graphical litera- 
ture a theorem on digraphs which contain a 2-basis. 





SOLUTIONS AND CLIQUES 

In order to develop an algorithm for finding all the solutions of a finite digraph, we 
require some additional concepts. Let D be a digraph, P be its set of points, and L its set 
of lines. If AC P, then the subgraph <A> generated by A has A as its set of points and 
contains all lines of D joining two points of A. Let Q be the collection of all ordered pairs 
of distinct points. Then the complement D of D is the digraph which also has P as its set of 
points, but whose set L of lines is given by Q- L. A subgraph of D is complete if there is a 
(directed) line of D joining every ordered pair of distinct points in it. A clique of a digraph is 
a maximal complete subgraph of three or more points. In [1] an algorithm for detecting all the 
cliques in any finite digraph is presented. We shall not reproduce this procedure here, but 
shall assume that it is available. Our next objective is to show the relationship between 
solutions and cliques. 

Let D be any digraph whose solutions (if any) we wish to find. We form the ordinary 
graph E from the complement D of D by deleting from D all those lines ab for which there 
does not occur the converse line ba. Note that D determines E uniquely, but not conversely. 
We now obtain theorems giving elementary criteria for internal stability and external stability. 


These will be combined to formulate a matrix algorithm for finding all the solutions of a finite 
digraph D. 


THEOREM 1. § is internally stable in D if and only if <S> is a complete subgraph 
of E. 

Proof: If S generates a complete subgraph of E, then any two points a, beS are 
adjacent in E. Hence both directed lines ab and ba are in D and § is internally stable in D. 


Conversely, if a, b€S, an internally stable set in D, then both ab and ba are in D. 
Thus <S> is a complete subgraph of E. 


The following observation now follows at once. 


COROLLARY 1. S is maximal internally stable in D if and only if S generates a 
maximal complete subgraph of E. 

But by definition, every solution of D is maximal internally stable in D since it is 
externally stable in D. Hence we have the additional result: 
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COROLLARY 2. Every solution of D generates a maximal complete subgraph of E. 


Since a clique has at least three points, Corollary 2 asserts that if S is a solution of D, 
then exactly one of the following possibilities hold in E: 

(1) S is an isolated point. 

(2) S consists of two adjacent points which do not lie in a triangle (3-cycle) of E. 

(3) <S> is a clique. 


Remark: This theorem and its corollaries are true for infinite digraphs as well, as the 
proof does not employ finiteness. 


We have just seen by the above theorem that all maximal complete subgraphs of E, or 
equivalently of D, are possible solutions of the finite digraph D. However, it is very easy to 
see that not every maximal complete subgraph of D, is necessarily a solution of D. An 
example is supplied by the graph D of Figure 1. Here we see from the graph E of Figure 1 
that the three possible solutions of D are {b}, {a, c}, or{a, d}. Referring to D, we find that 
{b} is not a solution since the line bd is not present. Among the remaining two possible solu- 
tions, {a, c} is a solution but {a, d} is not. 


d d 
d 
D: D: a b E: a eb 
a b 
c 
c c 
Figure 1 


Therefore, we need to add an additional condition (provided by Theorem 2) which a 
possible solution, i.e., the set of points of a maximal complete subgraph of E, must satisfy 
in order to be a solution of D. Those possible solutions which do not check this condition will 
be called extraneous. 

Let by, bo, beng b, be the points of D. Let A= (a;;) be the matrix for which aij = } 








if the line byb; is in D or if i=j, and aij = 0 otherwise; A is called the adjacency matrix 
of D. Let U=(1, 1, ... 1) be a vector with n components, and let R; be the ith row vector 


of A. By the addition of these row vectors, we will always mean boolean addition, 1+ 1 = 1. 


THEOREM 2. For notational convenience, let S={b,, bo, ..., bj}. Then S is 
externally stable in D if, and only if, 


Proof: Obvious. 
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Since the method in [1] for locating all the cliques in E (or equivalently in D) makes 
essential use of the fact that a clique has at least three points, we require the following: 
(1) The single point b; is a solution of D if and only if R, = U. 
Thus the singleton solutions of D can be read off from matrix A. 
(2) Two points b; and b; constitute a solution of D if and only if they are adjacent in 
E and R, + R; =U. 
It follows that b; and b; are not cocliqual in E, i.e., they do not lie in any clique. 
This is also implied by Corollary 2 of Theorem 1. 
(3) Finally, the set of points { by, Do, «+. b, where k> 2 isa solution of D if and 
only if this set is a clique in D and hence in E, and the following equation holds: 

















Ri, + Ro +... +R =U. 


Algorithm: 
For any given digraph D represented by adjacency matrix A = (a;;) whose ith row is 
R;; all the solutions (if any) may be found as follows: 
(1) Each row Ry, which equals U provides the singleton solution { b,}. ‘ 
(2) Any two points b, and b; of D adjacent in E form a solution { b,, bi} if R,+R;=U. 
(3) Find all the cliques of E using the straightforward method provided in [1]. Each 


such clique C is a solution if the equation 2 R, = U holds, where the sum is taken over the 
points of C. 


SUFFICIENT CONDITIONS FOR SOLVABILITY 

In the preceding section we have developed an algorithm for determining all the solu- 
tions of a finite digraph. This leaves unsolved the problem of characterizing solvable digraphs. 
There are, however, several sufficient conditions (Theorems 3, 4, and 5 below) each of which 
assures the solvability of a digraph. To state these, we require the following concepts. 

A (directed) cycle is a collection of k distinct points by, ree by, together with the lines 


by bo, evs Dye by, b, by - A semicycle also contains k distinct points and k directed lines, but 
here the lines are in the form: either of by bo or bob, -o-, either of b, by or b, 5. Thus 





every cycle is a semicycle, but not conversely. In [4], a semicycle is referred to as an 
“unoriented cycle." A cycle or semicycle is odd or even in accordance with its number of 
lines. A digraph is acyclic if it contains no directed cycles. It is strictly acyclic if it is 


acyclic and there are no infinite regressions [5]. The following three theorems are already 
known, 





THEOREM 3 (von Neumann and Morgenstern [3]). If D is acyclic, then D is solvable; 
further, D has a unique solution. 


THEOREM 4 (cf. [4]). If D has no odd semicycles, then D is solvable. 
THEOREM 5 (cf. [5]). If D has no odd cycles, then D is solvable. 


Clearly, Theorem 5 implies Theorem 4 and (excluding the uniqueness) Theorem 3 when 
D is finite. The relationship of the theorems in the infinite case is more subtle. If acyclicity 
is replaced by strict acyclicity, Theorem 3 remains true [3]. Theorem 4 remains true without 
modification [4]. If the additional hypothesis of local finiteness is made, Theorem 5 remains 
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true [5]. In [6] certain conditions are established for solvability of a digraph to be preserved 
under the operations of forming subgraphs and supergraphs. 

We now mention a few additional observations regarding the solvability of a digraph. 
The outdegree of a point b of D is the number of points adjacent from b; the indegree is 
defined similarly. A transmitter of a digraph is a point of indegree 0. 

(1) Every transmitter of a digraph lies in every solution. 

We note that this does not assert that every digraph with a transmitter necessarily has 
a solution; in fact this is not true. However, if S is a solution of D and b is a transmitter, 
then b €S. This is obvious since if b is not in § then there would exist a line of the form 
che D, which contradicts the hypothesis that b is a transmitter. 

A complete graph is one in which any two distinct points are adjacent. A complete 
oriented graph, or more briefly a tournament, is obtained from a complete graph by assigning 
a unique direction or orientation to each of its (unoriented) lines. With these definitions we 
obtain the following immediate consequence of the preceding observation. 

(2) A tournament is solvable if and only if it has a transmitter. 














Proof: The sufficiency is obvious since a transmitter in a tournament constitutes a 
singleton solution. To prove the necessity, we note that any solution of a complete oriented 
graph cannot contain more than one point; for between any two points of a tournament there is 
a directed line. But if a solution has exactly one point, then it must be a transmitter. 

Let T, (b) be the set of all points adjacent from b in D. We also call T,(b) the out-set 


of b. 
(3) If b is a transmitter of D, then D is solvable if and only if D-b-T ,(b) is solvable. 

This result, although trivial, may actually be useful. To solve a digraph which has 
transmitters, the first step (which serves to cut down the size of the digraph under considera- 
tion) is to delete all the transmitters and their out-sets. As an immediate generalization, 
we have: 

(4) If a point b of D has the property that it lies in every solution of D, then D is solvable if 
and only if D-b-T , (b) is solvable. 

A cycle of a digraph is maximal if its points do not lie in any longer cycle. As further, 
although more specialized, observations we have the following. 

(5) If D is a maximal odd cycle with exactly one diagonal line be, then D has a unique solu- 
tion containing b. — 

(6) Any digraph which consists of exactly one odd cycle has no solution. Any digraph which 
consists of exactly one even cycle has two solutions. 























THE CONCEPT OF r-BASIS 

A point b is accessible from c if there is a (directed) path from c to b. Let d(c, b) 
denote the distance from c to b, i.e., the length of any shortest path from c to b. A point 
basis is a minimal collection of points of the digraph G from which all points are accessible. 
The existence of a point basis for any finite digraph is shown in Kénig [2], and for infinite 
"ancestrally finite" digraphs in [5]. Arbitrary infinite digraphs need not have a point basis. 
We now show that there is a more general concept which includes both a point basis and a 
solution of a digraph as special cases. 


A subset S of the set P of points of a digraph G is an r-basis if the following two 
conditions hold: 
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(1) For any two distinct points of S, the distance from each to the other is greater than r. 
(2) For any point c € P - S, there exists be S such that d(b, c) =r. 

Conditions (1) and (2) may be termed internal and external r-stability, respectively. 

In these terms, a solution of a digraph is a 1-basis. Also, if n is the number of points of G, 
then any point basis of G is an (n-1)-basis. 

There already exists a theorem! which can be worded in terms of an r-basis for r = 2, 
namely: Every tournament has a singleton 2-basis. Translated into the language of pecking 
sequences of hens, this theorem asserts that in every community of hens there exists a hen 
such that every other hen is either pecked by her or is pecked by one which is pecked by her. 

We prove some theorems on the existence and uniqueness of r-bases in a digraph G. 














Given a digraph G, let an associated digraph cg) be formed using the same set of 
points as G but introducing a line xy in gt ) if and only if ye pb (x, G), 1S isr,y# x. The 


notation ye pi (x, G) means, as in [5], that there exists a directed path of length i in G from 
x to y. 


THEOREM 6. If § is a solution (or 1-basis) of go) then § is an r-basis of G, and 
conversely. 


Proof: Let S be a solution of cg), If x and y are points of S, then, by internal 
stability, there exists no line of g() joining them. Hence there exists no directed path of 
length less than or equal to r between them in G, so that S is internally r-stable in G. If y 
is a point not in S, then there exists a point x in S such that xy is in g) by external 
stability. Then y € D! (x, G), 15 i<r, so that S is externally r-stable in G. The converse 


is obvious. 


LEMMA, If G is acyclic, so is g), 





Proof: If there were a cycle of length k+1 in Gg), Say X, Xj Xp --- ee then 


i i i 
x, €D 1 (x), G), X,€D 2 (x, G),..-, x, €D ee G), where 15 i; =r for every j =1, 


k+1 
2, ..., k+l. Hence there exists a cycle of length ij in G, contradicting the hypothesis. 
j=1 


THEOREM 7. If G is finite and acyclic, then there exists an r-basis which is unique. 

Proof: This follows immediately from the preceding assertions of this section and the 
existence and uniqueness of a solution in a finite acyclic digraph established by von Neumann 
and Morgenstern in [3], which is Theorem 3 above. 


LEMMA. If G is strictly acyclic, so is cg), 
Proof: Strict acyclicity is equivalent to acyclicity together with regressive finiteness 
(cf. [5]). It is obvious that if G is regressively finite, so is G 54 





lSee for example Kemeny, J. G., J. L. Snell, and G. L. Thompson, Introduction to Finite 
Mathematics, Englewood Cliffs, 1957; pp. 317-318. 
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THEOREM 8. If G is infinite and strictly acyclic then there exists an r-basis which 
is unique. 
Proof: This again follows, as in Theorem 7, from the infinite version of Theorem 3, 


LEMMA. If G is locally finite, so is G), 
Proof: Immediate. 


THEOREM 9. If G is infinite, locally finite, and acyclic (but not necessarily strictly 
acyclic), then there exists an r-basis. 

Proof: This follows from the preceding assertions together with the corollary to 
Theorem 11 of [5]. 


Uniqueness does not hold in this latter case as simple examples show. 

By virtue of Theorem 6, the algorithm for finding the solutions or 1-bases for a given 
finite digraph G can be used to find the r-bases by finding the solutions of the associated 
digraph cg), The adjacency matrix of go ) can easily be calculated from A, the matrix of 
G; in fact it is obtained by forming At using boolean addition. 

That Theorem 6 goes beyond the acyclic case is shown by the following example: 


b b 


Gs: g!2) : 


Figure 2 


In Figure 2, {a} is a 2-basis of G and a solution of c(2). phat not every digraph has an r-basis 
can be seen from the following simple example possessing no 2-basis. 
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Remark. The theorems of [4, 5, 6] concerning the absence of cycles of odd length 
cannot be extended to similar theorems on r-bases because of the following proposition: If 
G has any cycle of length greater than two, then gf ). r> 1, has an odd cycle. This follows 
from the fact that if r = 2 and m = 3 then there exists a partition of the integer m with an 
odd number of summands, m = =m,, where all m,; =r. For we may write m=qr+s, 
0=s<pr. If q is odd and s = 0, there is nothing to prove. If q is even and s > 0, there is 
again nothing to prove. If q is even and s= 0, we may write m = (q-1)r + (r-1) +1. If q is 
odd and s > 0, we may write m = (q-1)r + s + (r-1) +1. 

Problem 1 of the introduction may now be generalized as follows: Characterize those 


irreflexive relations for which an r-basis exists, and those for which a unique r-basis exists, 
for arbitrary r2 1. 








SUMMARY 


An irreflexive binary relation, such as domination of imputations in n-person games 
as developed by von Neumann and Morgenstern, is represented by a digraph (directed graph). 
A set § of points of the digraph is a solution if (1) no point of S is dominated by another point 
of S and (2) every point not in S is dominated by some point of S. A matrix algorithm is 
developed for finding all the solutions, if any, of a finite irreflexive relation. First all cliques 
(maximal complete subgraphs) of the complementary relation are found. For such a clique to 
be a solution, it must then satisfy a further matrix condition, which is easily checked. The 
concept of solution is generalized as follows. An r-basis of a digraph is a set S of points 
such that: 

(1) no two points of S are joined by a directed path from one to the other of length less than or 
equal to r; (2) every point not in S is the end of a directed path of length less than or equal to 

r with initial point in S. It is proved that a digraph has an r-basis if it is (1) finite and acyclic, 
or (2) infinite and strictly acyclic, or (3) infinite, locally finite, and acyclic but not necessarily 
strictly acyclic. For a finite digraph, the matrix algorithm for solutions can be extended to 
obtain the r-bases. 


NOTE: The authors have just learned that the recent book C. Berge, Theorie des graphes et 


ses applications, Dunod, Paris, 1958, contains an independent and different treatment for 
finding the solutions of a directed graph. 


REFERENCES 


[1] Harary, F. and I. C. Ross, "A Procedure for Clique Detection Using the Group Matrix," 
Sociometry 20:205-215 (1957). 


K6nig, D. Theorie der endlichen und unendlichen Graphen, reprinted Chelsea Publishing 
Company 1950. 





von Neumann, J. and O. Morgenstern. Theory of Games and Economic Behavior, Princeton 
University Press, 1944; 2nd edition, 1947. 





Richardson, M., "On Weakly Ordered Systems," Bull. Amer. Math. Soc. 52:113-116 (1946). 





, "Solutions of Irreflexive Relations,"" Annals of Mathematics, 58:573-590 (1953). 





, "Relativization and Extension of Solutions of Irreflexive Relations," Pacific 
Journal of Mathematics, 50:551-584 (1955). 





* * 








x 








NOTES 


PROBLEMS 

It has been suggested that this journal might serve as a medium for publishing problem 
in the area of Logistics — with the idea that other persons might be interested in those proble1 
and might submit comments. Readers are invited to submit brief statements on applied and 
theoretical problems in Logistics. Address letters to Managing Editor, Naval Research Logis 
tics Quarterly, Office of Naval Research, Washington 25, D. C. 


ERRATA 


The title of the paper "The New Military Standard for Inspection by Variables" which 
appeared in the June 1959, Volume 6, Number 2 issue of this journal should have been "The 
New Military Standard 414 for Acceptance Inspection by Variables." 
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NEWS AND MEMORANDA 


Readers are invited to submit to the Managing Editor items of 
general interest in the field of logistics 


The second West Coast ORSA-TIMS meeting will be held at the U. S. Naval Postgraduate 
School, Monterey, California on April 7 and 8, 1960. The meeting is sponsored jointly by the 
Northern and Southern California Chapters of TIMS and the Western Section of ORSA. There 
will be both invited and contributed sessions of theoretical and applied character. Abstracts of 
contributed papers must reach the Program Chairman by February 1, 1960. The committee 
chairmen for this meeting are: 


General Chairman - Prof. D. Teichroew, 
Stanford University 


Arrangements - Prof. T. Oberbeck, 
Chairman U. S. Naval Postgraduate School 


Program Chairman - E. Koenigsberg 
Touche, Niven, Bailey & Smart 

100 Bush Street 

San Francisco 4, California 
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MANAGEMENT TRAINING -- Cases and Principles, by William J. McLarney. 
Richard D. Irwin, Inc., Homewood, I1l., 1959. xvii + 521 pp., $8.10. 

This is the third edition of a manual recommended by the publishers for use as a text 
for students of management or personnel in supervisory training courses. The main unique 
feature of the book is a series of 286 cases. These cases are usually quite brief, averaging 
about 250 words each. They are arranged by topic, and follow generalized discussion of 
principles bearing on the topic. Each case is followed from one to eight questions about the 
application of the principles to the case. Some of the questions are quite general, such as 
"How do you view the situation?" or "What can you do to correct the situation?" or 'How 
would you handle the situation?" 

The four major parts of the book are: 

--Management's Job (including sections on organization policies and procedures, chain 
of command and cooperation). 

--Job Management (including planning, methods, standards of performance and place- 
ment), 








--Developing the Work Team (including inducting, evaluating and transferring 
personnel). 

--Maintaining the Work Team (including supervisory leadership, discipline, grievances 
and morale). 

Each section is followed by a series of questions designed to provoke interest in follow- 
ing up on the job. 

A generalized set of steps for problem solving is a main feature. The steps recom- 
mended by the author are: 

1, -‘Clearly define the problem or problems. 

2. Gather the information. 

3. Interpret the information. 
4. Develop solutions. 
5. 
6. 








Select the best practical solution. 
Put the solution into effect. 

7. Evaluate the effectiveness of the solution. 

These steps are considered the proper procedure to be applied to the case problems. 
There may be some lack of emphasis on testing of alternative solutions in the foregoing seven- 
step method. 

Mr. McLarney places considerable emphasis on the necessity for supervisors managing 
the work as well as managing people. In this sense the book is well-balanced. His section on 
induction contains some excellent hints to supervisors; his section on evaluation of perform- 
ance is weakened by insistence that it is almost always good to tell employees almost 
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everything about their performance. This is a highly debatable point in terms of proper 
motivation and employee-supervisory relationships in many situations. 

This book is a textbook. It is not for casual reading at one sitting. In view of the 
generalized nature of many of the cases and questions, it is probably more useful as a text for 
class discussion, than for self-study. 

(Reviewed by C. E. Auvil) 


LOGISTICS IN THE NATIONAL DEFENSE, by Rear Admiral Henry E. Eccles, U. S. 
Navy (Retired). The Stackpole Company, Harrisburg, Pennsylvania, 1959. 347 pp. $5.00. 

Logistics in the National Defense deals not only with the practice of logistics at the 
higher levels of military command and civilian administration, but with the theory of logistics. 

Admiral Eccles points out that we are in the midst of an electronic-nuclear revolution 
during which economics is the limiting factor in the development of our military defense sys- 
tem. By this backdoor, every taxpayer becomes vitally interested in the capabilities of our 
military Services and our overall strength and national security. 

Admiral Eccles believes that the actions and decisions of civilian control and military 
command, at all levels, are bound to be based upon a blend of strategical, logistical and tactical 
considerations. He develops by example how strategy is shaped by economic and logistic 
capabilities. 

The author then points out how, in the Pentagon, there has been a tendency to subordinate 
the basic principles of the military art to the terminology and practices of a business world, in 
which the basic criteria are quite different from the criteria of military excellence or of suc- 
cess in combat operations. 

With this background, Admiral Eccles expresses the belief that any indifference to 
logistic organization or planning, or any incompetence in the planning and supervision of 
logistical activities on the part of senior line officers of our three Services is an open invitation 
to the further erosion, or perhaps even the destruction of military command authority in broad 
areas of operational logistics. 

The various components of modern conflict—strategical, logistical and tactical—are 
examined. Each component is related to the other and to the whole spectrum of modern war. 
How strategy is shaped by economic and logistic capabilities is developed by examples. 

The problem, constantly encountered, of logistic activities tending to snowball or to 
grow out of all proportions to the tactical forces which they support is dealt with fully, fre- 
quently, and with marked perception. 

Admiral Eccles believes that sound logistics form the foundation for the development of 
strategic flexibility and mobility. Hence, military commanders must know and control their 
logistics. It is pointed out, with examples, where the strong tendency to control operational 
practices from the higher reaches of the Pentagon creates a rigidity which defeats the purposes 
of control, and that similarly, overly large staffs tend to obscure the commanders personal 
responsibility in the process of integrated planning. 

Whether or not you agree with Admiral Eccles that disarmament in the nuclear age is 
even a more remote possibility than it has been in the past, that Suez was a logistic defeat for 
the British, that inefficient personnel is the greatest single source of the logistic snowball, that 
movement control should operate as an individual staff agency or subordinate command, or 
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that it will be largely in the field of logistics that our readiness for future conflict will be 
determined, you will find reasoned support for each of these statements and many others 
equally challenging. 

Logistics in the National Defense is no funny picture book with underdigested state- 
ments and oversimplified conclusions. It is not for reading on the run. But it is an important 
book, a tremendously stimulating book, a reasoned book. It is a landmark for logisticians, 
marking the most comprehensive and worthwhile advance in the study of the theory of logistics, 
since Colonel Cyrus Thorpe USMC wrote his "Pure Logistics" in 1917. It is professional 
military education at its timely best, and yet profitable reading for the layman as is so ably 
pointed out by Dr. Henry G. Wriston, President of the American Assembly, Columbia University, 
in the Foreword. 

(Reviewed by G. C. Dyer) 


THE QUESTION OF NATIONAL DEFENSE, by Oskar Morgenstern. Random House, 
New York, 1959. 296 pp. 

The distinguished Professor of Political Economy at Princeton and co-author of the 
Theory of Games, after many years of association with the scientific aspects of the problem, 
has written a frank, provocative, and at times bitterly ironical analysis of national defense. 

After a fine prologue and a chapter on the present threat, Professor Morgenstern 
continues with estimates of possible reactions and a discussion of provocation. -The nature 
and importance of stability in relation to threat, and a brilliant analysis of the theory of 
deterrence are among the many fresh ideas expressed. 

He explains and stresses the distinction between counter-attack and retaliation, 
particularly as related to vulnerability. 

After commenting on the possibilities of accidential action, he moves on to an extensive 
discussion of vulnerability which brings out the merits of a mixed system retaliatory force 
and the fallacy of major dependence on “hardened” missile bases. This leads to a complete 
statement of his concept of the Oceanic System previously published in September 1959 
Fortune. 

The key to the great value of the Oceanic System lies in the thought that the capability 
of selective, carefully timed nuclear retaliation is a better deterrent than the threat of 
massive simultaneous use of all nuclear weapons. 

The tremendous advantages in concealment, mobility, and unpredictability of the 
Polaris submarine concept are now, of course, well recognized. The suitability of his nuclear 
seaplane concept, however, is debatable. 

As a result of his extensive thought on game theory, Professor Morgenstern makes the 
further vital, but heretofore little understood, point that it is to our advantage for Russia to 
to have an invulnerable force. This is clearly stated in two basic principles and a further 
comment. 

“In order to preserve a nuclear stalemate under conditions of nuclear 
plenty it is necessary for both sides to possess invulnerable retaliatory 
forces. (p. 74) 
“In view of modern technology of speedy weapons delivery from any 

poi point on earth to any other, it is in the interest of the United States for 
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Russia to have an invulnerable retailiatory force and vice versa. 

(p. 75-76) 

“There should be no misunderstanding of what it would mean for Russia 
to have an invulnerable retailiatory force. We would lose nothing there- 
by-absolutely nothing except the ability to stage a successful surprise 
attack against the enemy.” (p. 77) 

There are excellent discussions of the limitation of war, and the economic burden which 
can be sustained. The comments on “military worth” are particularly valuable. 

Chapter V is of such vital importance that it deserves special comment. The chapter is 
devoted to the difference between wars of attrition and those of annihilation; an estimate of the 
casualties to be expected and a plea for the construction of shelters against “fall out.” His 
comment on the inadequate concepts and inept administration of the U.S. Civil Defense Program 
is scathing and spares neither the President, the Congress, nor the administrators! 

This reviewer would hope that the author dwell on the further and paradoxical implica- 
tions of this very important chapter. The first is: What will be the effect of such shocking 
casulaties and destruction on the political and social concepts of the survivors? The second: 
What type of civilian-military control will be required to manage the remnants of our society 
and to direct the enormous and crude logistical operations of recovery? The third is: Since 
these remnants cannot re-establish the civilized ideas and concepts which we supposedly 
fought to preserve except by the most rigid discipline and dedicated purpose, what would be 
the effect on our military position today were such sense of responsibility and purpose achieved? 

It would seem likely that if our people can somehow acquire this sense of purpose and 
discipline now, thermonuclear war can be avoided. Certainly such purpose and discipline would 
enormously strengthen our military economic and political position. The extreme measures of 
large program of shelter construction and placing part of the economy underground would then 
not be needed. 

Morgenstern touches on these matters on page 129 where he says: 

“Recovery is therefore something much more elemental. It aims at 
continuation of life, on a simple and primitive basis. From mere sur- 
vival, in view of perhaps a hundred million people killed within a few 
days, to the construction of a society which adheres to the lofty ideas 
we cherish today is a long way. And the survivors may develop ideals 
so different from ours that we may not recognize the country— should 
we live to make the comparison, or, living, should we care to make it 
even if we could. 

“Still, we now believe that some values of our civilization should be 
preserved. Therefore we must make this possible. At present we do 
only a partial job poorly.” 

The latter part of his book discusses the security process, information, intelligence, 
politics and diplomacy~—all as they involve defense. This part includes an expression of 


appreciation for the problems of high military command, and a plea for a more thorough and 
scientific study of politics. 


Finally, toward the close there is a perceptive and important comment on negotiations 
with the enemy. This should be required reading for all devotees of the push button concept! 
This challenging book ends with the pessimistic comment: 


























“As it is, the probability of a large thermonuclear war occurring appears 
to be significantly larger than the probability of its not occurring. Will 
at least these probabilities be reversed ?” 

Very definitely Professor Morgenstern is not happy with many of our official policies and 
is particularly displeased with official attitudes. He comments on them forcefully and with good 
reason. His frequently expressed disdain for superficiality and cant, and his pleas for intellec- 
tual rigor are indeed refreshing. 

Recognizing the great knowledge and dedication which has gone into the writing of this 
important work, this reviewer suggests the need for a further appreciation. 

Specifically, while the great danger of thermonuclear war should not be minimized, the 
problem can be put into its correct perspective only when there has first been a careful study 
of the nature of the spectrum of human conflict, an analysis (not merely a statement) of national 
objectives and, most important of all, a thorough analysis of the nature of strategy itself. 

Until these are done, even such excellent books as this will be inadequate and perhaps 
only tentative approaches to the vital questions of our national defense. 


(Reviewed by H. E. Eccles) 
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1954, pp. 230-232. 





*In this index, no extensive attempt has been made to cover the many possibilities, such as: 
"falls under," "is related to,'"' ''cites a connection with," ''makes application to," etc. In 
some cases, papers were placed in several categories because of connections notimmediately 
apparent from their titles. Conversely, some titles are so clear as to preclude the necessity 
for certain subheadings (or even cross references) or multiple listings. 
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Decision Making—Continued 
"Decision Making in the Face of Uncertainty—II"' - R. Bellman, Vol. 1, No. 4, December 
1954, pp. 327-332. 
"Organized Decision Making" - H. D. Mills, Vol. 2, No. 3, September 1955, pp. 137-143. 
"Decision Making and New Mathematics" - G. L. Thompson, Vol. 3, No. 3, September 
1956, pp. 141-149. 
"A Formulation of the Decision Problem for a Class of Systems" - J. K. Hale and R. L. 
Reed, Vol. 3, No. 4, December 1956, pp. 259-277. 
"Linear Model for Evaluating Complex Systems" - R. M. Thrall, C. H. Coombs, and 
W. Caldwell, Vol. 5, No. 4, December 1958, pp. 347-361. 
"Allocation of a Resource to Alternative Probabilistic Demands: Transport-Equipment 
Pool Assignments" - R. Saposnik and V. L. Smith, Vol. 6, No. 3, September 1959, 
pp. 193-207. 
Demand, Analysis of (see Inventory and Supply, Demand and Requirements Analysis) 
Demand Distribution (see Inventory and Supply, Demand and Requirements Analysis) 
Demand, Low (see Inventory and Supply, Demand and Requirements Analysis) 
Dental Materials (see Medical) 
Development Funds (see Economic Models and Problems, Managerial Economics) 
Discrete Programming (see Programming, Linear) 
Distribution, Demand (see Inventory and Supply, Demand and Requirements Analysis) 
Distribution of Inventories (see Inventory and Supply, Distribution) 
Dual Algorithm (see Programming, Linear) 
Dual Simplex Method (see Programming, Linear) 
Duels (see Games Models) 
Dynamic Inventory Problems (see Inventory and Supply, Demand and Requirements Analysis) 
Dynamic Programming (see Programming, Dynamic) 


Economic Models and Problems, General 
"On the Graphs and Asymptotic Forms of Finite Boolean Relation Matrices and Stochastic 
Matrices" - D. Rosenblatt, Vol. 4, No. 2, June 1957, pp. 151-167. 
"Dynamic Problems in the Theory of the Firm" - H. M. Wagner and T. M. Whitin, Vol. 5, 
No. 1, March 1958, pp. 53-74; also see Editor's Note, Vol. 5, No. 2, June 1958, 
p. 189. 
Economic Models and Problems, Industrial Mobilization and Planning 
"On the Meaning and Use of a Capacity Concept" - E. Zabel, Vol. 2, No. 4, December 1955, 
pp. 237-249. 
"Measures of Industry Capacity'' - E. Zabel, Vol. 3, No. 4, December 1956, pp. 229-244. 
"Damage Assessment Systems and Their Relationship to Post-Nuclear-Attack Damage 
and Recovery" - F. A. Ramsey, Jr., Vol. 5, No. 3, September 1958, pp. 199-219. 
Economic Models and Problems, Input-Output Systems 
"An Elementary Treatment of an Input-Output System" - Y. K. Wong, Vol. 1, No. 4, Dec. 
1954, pp. 321-326. 
"Probabilistic Errors in the Leontief System" - R. E. Quandt, Vol. 5, No. 2, June 1958, 
pp. 155-170. 
"On the Solution of Probabilistic Leontief Systems" - R. E. Quandt, Vol. 6, No. 4, Decem- 
ber 1959, pp. 295-305. 
Economic Models and Problems, Managerial Economics 
"Linear Programming in Bid Evaluation" - E. D. Stanley, D. P. Honing, and L. Gainen, 
Vol. 1, No. 1, March 1954, pp. 48-54. 
"Contractor Performance" - D. C. MacKenzie, Vol. 2, No. 1 & 2, March & June 1955, 
pp. 17-23. 
"A Bound on the Use of Inefficient Indivisible Units"' - M. J. Beckmann and J. Laderman, 
Vol. 3, No. 4, December 1956, pp. 245-252. 
"Experiences with the Bid Evaluation Problem" (Abstract) - H. Bremer, W. Hall, and 
M. Paulsen, Vol. 4, No. 1, March 1957, p. 27. 
“Computing Procedures for Portfolio Selection" (Abstract) - H. M. Markowitz, Vol. 4, 
No. 1, March 1957, pp. 87-88. 
"A Transportation and Production Model" - G. W. Evans, II, Vol. 5, No. 2, June 1958, 
pp. 137-154. 
"The Allocation of Development Funds: An Analytic Approach" - S. C. Daubin, Vol. 5, No. 3, 
September 1958, pp. 263-276; see also Errata Sheet, Vol.6, No. 2, June 1959, pp. 185-186. 
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Economic Models and Problems, Process Analysis 
"Linear Programming Planning of Refinery Operations" - L. K. Cheney, Vol. 4, No. 1, 
March 1957, pp. 9-16. 
Economic Models and Problems, Utility 
"Priorities in the Naval Supply Systems" - W. M. Young, Vol. 1, No. 1, March 1954, 
pp. 16-24. 
"Optimal Weapon Systems" - J. Marschak and M. R. Mickey, Vol. 1, No. 2, June 1954, 
pp. 116-140. 
"A Plan to Allocate and Procure Electronic Sets by the Use of Linear Programming 
Techniques and Analytical Methods of Assigning Values to Qualitative Factors" - 
J. W. Smith, Vol. 3, No. 3, September 1956, pp. 151-162. 
"Procurement and Allocation of Naval Electronic Equipments" - G. Suzuki, Vol. 4, No. 1 
March 1957, pp. 1-7. 
"A Selected, Descriptive Bibliography of References on Priority Systems and Related, 
Nonprice Allocators"' - W. G. Mellon, Vol. 5, No. 1, March 1958, pp. 17-27. 
"A Method of Estimating Spare-Part Essentiality" - H. W. Karr, Vol. 5, No. 1, March 
1958, pp. 29-42. 
"The Coefficients in an Allocation Problem" - R. J. Aumann and J. B. Kruskal, Vol. 5, 
No. 2, June 1958, pp. 111-123. 
"Notes" - R. E. McShane and H. Solomon, Vol. 5, No. 4, December 1958, pp. 363-367. 
"Assigning Quantitative Values to Qualitative Factors in the Naval Electronics Problem" - 
R. J. Aumann and J. B. Kruskal, Vol. 6, No. 1, March 1959, pp. 1-16. 
"Notes" - R. G. Davis, Vol. 6, No. 2, June 1959, pp. 183-186. 
"Notes" - J. B. Kruskal, Vol. 6, No. 3, September 1959, pp. 261-262. 
Electronic Equipments, Allocation and Procurement of (see Economic Models and Problems, 
Utility ) 
Employment Scheduling (see Scheduling) 
Essentiality (see Economic Models and Problems, Utility; also Inventory and Supply, Demand 
and Requirements Analysis) 
Excess Stock (see Inventory and Supply, Stocking Policies) 


Firm, Theory of (see Economic Modeis and Problems, General) 


Game Theory 

"A Numerical Method to Determine Optimum Strategy" - J. von Neumann, Vol. 1, No. 2, 
June 1954, pp. 109-115. 

"Fictitious Play for Continuous Games" - J. M. Danskin, Vol. 1, No. 4, December 1954, 
pp. 313-320. 

"The Concept of a 'Mission'' - M. Verhulst, Vol. 3, No. 1 & 2, March & June 1956, 

. 45-57. 

"Decision Making and New Mathematics" - G. L. Thompson, Vol. 3, No. 3, September 
1956, pp. 141-149. 

"The Basic Theorems of Real Linear Equations, Inequalities, Linear Programming, and 
Game Theory" - D. Gale, Vol. 3, No. 3, September 1956, pp. 193-200. 

"Some Methods of Computational Attack on Programming Problems Other than the Simplex 
Method" (Abstract) - C. Thompkins, Vol. 4, No. 1, March 1957, pp. 95-96. 

"Equilibrium Points in Games with Vector Payoffs'' - L. S. Shapley (Introductory Note by 
F. D. Rigby), Vol. 6, No. 1, March 1959, pp. 57-61. . 

"A Terminating Iterative Algorithm for Solving Certain Games and Related Sets by Linear 
Equations" - R. B. Braithwait, Vol. 6, No. 1, March 1959, pp. 63-74. 

Games Models 
"Some Blotto Games" - D. W. Blackett, Vol. 1, No. 1, March 1954, pp. 55-60. 
"On Multi-Component Attrition Games" - D. Blackwell, Vol. 1, No. 3, September 1954, 


pp. 210-216. 

"Decision Making in the Face of Uncertainty—II" - R. Bellman, Vol. 1, No. 4, Dec. 1954, 
pp. 327-332. 

"The Problem of Aiming and Evasion" - R. Isaacs, Vol, 2, No. 1 & 2, March & June 1955, 
pp. 47-67. 


"Analytic Study of War Games" - H. A. Antosiewicz, Vol. 2, No. 3, September 1955, 
pp. 181-208. 
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Games Models—Continued 
"Attrition Games" - J. R. Isbell and W. H. Marlow, Vol. 3, No. 1 & 2, March & June 1956, 
pp. 71-94. 
"A Tank Duel with Game-Theoretic Implications" - L. E. Zachrisson, Vol. 4, No. 2, June 
1957, pp. 131-138. 
"An Application of Game Theory to Special Weapons Evaluation" - J. K. Hale, and H. H. 
Wicke, Vol. 4, No. 4, December 1957, pp. 347-356. 
"Pure Strategy Solutions of Blotto Games" - D. W. Blackett, Vol. 5, No. 2, June 1958, 
pp. 107-109. 
"Games of Economic Survival" - M. Shubik and G. L. Thompson, Vol. 6, No. 2, June 1959, 
pp. 111-123. 
"The Watchdog and the Burglar - S. P. Thompson and A. J. Ziffer, Vol. 6, No. 2, June 
1959, pp. 165-172. 
Graph Theory (see Miscellaneous Mathematics) 


Hitchcock-Koopmans Problem (see Programming, Linear—The Transportation Problem) 
Hotspot (see Games Models) 


Hungarian Method (see Programming, Linear—The Assignment Problem) 


Industrial Capacity (see Economic Models and Problems, Industrial Mobilization and Planning) 
Industrial Mobilization and Planning (see Economic Models and Problems, Industrial Mobiliza- 
tion and Planning) 
Input-Output Systems (see Economic Models and Problems, Input-Output Systems) 
Inspection by Variables (see Reliability and Quality Control) 
Intelligence Systems (see Data Processing) 
Inventory and Supply, General 
"Cycling" - I. R. Savage, Vol. 3, No. 3, September 1956, pp. 163-175. 
"An Inventory Control Bibliography" - R. Lewis, F. Neeland, and M. Gourary, Vol. 3, 
No. 4, December 1956, pp. 295-303. 
"Renewal-Theoretic Analysis of a Two-Bin Inventory Control Policy" - D. P. Gaver, Jr., 
Vol. 6, No. 2, June 1959, pp. 141-163. 
Inventory and Supply, Demand and Requirements Analysis 
"Statistical Methods for Determining Requirements of Dental Materials" - J. A. English 
and E. A. Jerome, Vol. 1, No. 3, September 1954, pp. 191-199. 
"Science and Logistics" - R. E. McShane, Vol. 2, No. 1 & 2, March & June 1955, pp. 1-7. 
"A Method for Calculating Optimal Inventory Levels and Delivery Time" - T. M. Whitin 
and J. W. T. Youngs, Vol. 2, No. 3, September 1955, pp. 157-173. 
"Warehousing and Distribution of a Seasonal Product" - W. S. Jewel, Vol. 4, No. 1, March 
1957, pp. 29-34. 
"A Method of Estimating Spare-Part Essentiality" - H. W. Karr, Vol. 5, No. 1, March 
1958, pp. 29-42. 
“Airline Demand: An Analysis of Some Frequency Distributions" - M. J. Beckmann and 
F. Bobkoski, Vol. 5, No. 1, March 1958, pp. 43-51. 
"Loss Resulting from the Use of Incorrect Data in Computing an Optimal Inventory 
Policy" - J. Levy, Vol. 5, No. 1, March 1958, pp. 75-81. 
"Dynamic Problems in the Theory of the Firm" - H. M. Wagner and T. M. Whitin, Vol. 5, 
No. 1, March 1958, pp. 53-74. (See Notes, Vol. 5, No. 2, June 1958, pp. 189) 
"A Transportation and Production Model" - G. W. Evans, Vol. 5, No. 2, June 1958, 
pp. 137-154. 
"Characteristics of Usage of Supply Items Aboard Naval Ships and the Significance to 
Supply Management" - J. H. Garrett, Vol. 5, No. 4, December 1958, pp. 287-306. 
"Further Notes on the Loss Resulting from the Use of Incorrect Data in Computing an 
Optimal Inventory Policy" - J. Levy, Vol. 6, No. 1, March 1959, pp. 25-31. 
"Renewal-Theoretic Analysis of a Two-Bin Inventory Control Policy" - D. P. Gaver, Jr., 
Vol. 6, No. 2, June 1959, pp. 141-163. 
"Allocation of a Resource to Alternative Probabilistic Demands: Transport-Equipment 
Pool Assignments" - R. Saposnik and V. L. Smith, Vol. 6, No. 3, September 1959, 
pp. 193-207. 
"A Method for Determining Supply Quantity for the Case of Poisson Distribution of 
Demand" - M. R. Mickey, Jr., Vol. 6, No. 4, December 1959, pp. 265-272. 
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Inventory and Supply, Distribution 
"Procurement and Allocation of Naval Electronic Equipments" - G. Suzuki, Vol. 4, No. 1, 
March 1957, pp. 1-7. 
"Warehousing and Distribution of a Seasonal Product" - W. S. Jewell, Vol. 4, No. 1, March 
1957, pp. 29-34. 
"Redistribution of Total Stock Over Several User Locations" - S. G. Allen, Vol. 5, No. 4, 
December 1958, pp. 337-345. 
Inventory and Supply, Management 
"Priorities in the Naval Supply System" - W. M. Young, Vol. 1, No. 1, March 1954, 
pp. 16-24. 
"On the Span of Central Direction" - T. M. Whitin, Vol. 1, No. 1, March 1954, pp. 25-35. 
"Navy Stratification and Fractionation for Improvement of Inventory Management" - 
F. L. Hetter, Vol. 1, No. 2, June 1954, pp. 75-78. 
"An Application of Statistical Techniques to Management of Overseas Supply Operations"' 
E. F. Hollingshead, Vol. 1, No. 2, June 1954, pp. 82-89. 
"Developing an Engineering Productivity Standard" - R. C. Horne, Vol. 1, No. 3, Septem- 
ber 1954, pp. 203-206. 
"Consistency Problems in the Military Supply System" - O. Morgenstern, Vol. 1, No. 4, 
December 1954, pp. 265-281. 
"The Theory of Inventory Management—A Review" - T. V. V. Atwater, Jr., Vol. 1, No. 4, 
December 1954, pp. 295-300. 
"The Costs of Alternative Air Base Stocking and Requisitioning Policies''- J. W. Petersen 
and M. A. Geisler, Vol. 2, No. 1 & 2, March & June 1955, pp. 69-81. 
"A Formula for Decisions on Retention or Disposal of Excess Stock" - J. R. Simpson, 
Vol. 2, No. 3, September 1955, pp. 145-155. 
"Comments on the Morgenstern Model" - J. C. Busby, Vol. 2, No. 4, December 1955, 
pp. 225-236. , 
"A Selected, Descriptive Bibliography of References on Priority Systems and Related, 
Nonprice Allocators" - W. G. Mellon, Vol. 5, No. 1, March 1958, pp. 17-27. 
"Some Principles for a Data-Processing System in Logistics"’- M. A. Geisler, Vol. 5, 
No. 2, June 1958, pp. 95-105. 
"Characteristics of Usage of Supply Items Aboard Naval Ships aud the Significance to 
Supply Management" - J. H. Garrett, Jr., Vol. 5, No. 4, December 1958, pp. 287-306. 
Inventory and Supply, Procurement 
"Design Change Impacts on Airframe Parts Inventories" - J. W. Petersen and W. A. 
Steger, Vol. 5, No. 3, September 1958, pp. 241-255. 
Inventory and Supply, Production and Sales 
"A Scientific Method for Establishing Reorder Points'' - M. J. Solomon, Vol. 1, No. 4, 
December 1954, pp. 289-294. 
"A Multi-Stage Inventory Model" - J. G. Bryan, G. P. Wadsworth, and T. M. Whitin, 
Vol. 2, No. 1 & 2, March & June 1955, pp. 25-37. 
"Mathematical Treatment of a Stockpiling Problem" - J. M. Danskin, Vol. 2, No. 1 & 2, 
March & June 1955, pp. 99-109. 
"Warehousing and Distribution of a Seasonal Product" - W. S. Jewell, Vol. 4, No. 1, 
March 1957, pp. 29-34. 
"Some Production Planning Problems" - M. Klein, Vol. 4, No. 4, December 1957, 
pp. 269-286. 
"An inventory Policy for Repair Parts" - M. J. Beckmann, Vol. 6, No. 3, September 1959, 
pp. 209-220. 
Inventory and Supply, Qualitative Aspects 
"Issue Priority: Last In First Out (LIFO) vs First In First Out (FIFO) as a Method of 
Issuing Items from Supply Storage" - J. A. Greenwood, Vol. 2, No. 4, December 
1955, pp. 251-268. 
Inventory and Supply, Replenishment Systems 
"Priorities in the Naval Supply System" - W. M. Young, Vol. 1, No. 1, March 1954, 
pp. 16-24. 
"The Application of Statistical Techniques to Management of Overseas Supply Operations” - 
E. F. Hollingshead, Vol. 1, No. 2, June 1954, pp. 82-89. 
"Mobile Logistics Support in the 'Passage to Freedom' Operation" - R. Meade, Jr., and 
C. A. Fischer, Vol. 1, No. 4, December 1954, pp. 258-264. 
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Inventory and Supply, Replenishment Systems—Continued 
"The Costs of Alternative Air Base Stocking and Requisitioning Policies" - J. W. Petersen 
and M. A. Geisler, Vol. 2, No. 1 & 2, March & June 1955, pp. 69-81. 
"Multinational Logistics in the Nuclear Age" - C. L. Henn, Jr., Vol. 4, No. 2, June 1957, 
pp. 117-129. 
"Improvise or Plan?" - A. H. Wilson and W. R. Finn, Vol. 4, No. 4, December 1957, 
pp. 263-267. 
Inventory and Supply, Stocking Policies 
"Notes on Applied Analytical Logistics in the Navy" - C. F. Lynch, Vol. 1, No. 2, June 
1954, pp. 90-102. 
"The Costs of Alternative Air Base Stocking and Requisitioning Policies" - J. W. Petersen 
and M. A. Geisler, Vol. 2, No. 1 & 2, March & June 1955, pp. 69-81. 
"Mathematical Treatment of a Stockpiling Problem" - J. M. Danskin, Vol. 2, No. 1 & 2, 
March & June 1955, pp. 99-109. 
"A Formula for Decisions on Retention or Disposal of Excess Stock" - J. R. Simpson, 
Vol. 2, No. 3, September 1955, pp. 145-155. 
"A Method for Calculating Optimal Inventory Levels and Delivery Time" - T. M. Whitin 
and J. W. T. Youngs, Vol. 2, No. 3, September 1955, pp. 157-173. 
"An Optimum Allowance List Model" - M. H, Gourary, Vol. 3, No. 3, September 1956, 
pp. 177-191. 
"A Simple Rule for the Consolidation of Allowance Lists" - M. H. Gourary, Vol. 5, No. 1, 
March 1958, pp. 1-15. 
"Loss Resulting from the Use of Incorrect Data in Computing an Optimal Inventory 
Policy" - J. Levy, Vol. 5, No. 1, March 1958, pp. 75-81. 
"Further Notes on the Loss Resulting from the Use of Incorrect Data in Computing an 
Optimal Inventory Policy" - J. Levy, Vol. 6, No. 1, March 1959, pp. 25-31. 
"A Method for Determining Supply Quantity for the Case of Poisson Distribution of 
Demand" - M. R. Mickey, Jr., Vol. 6, No. 4, December 1959, pp. 265-272. 
Inventory and Supply, Warehousing 
"Warehousing and Distribution of a Seasonal Product" - W. S. Jewell, Vol. 4, No. 1, 
March 1957, pp. 29-34. 
Issue Priority (see Inventory and Supply, Qualitative Aspects) 


Job Lot Scheduling (see Scheduling) 
Job Shop Scheduling (see Scheduling) 


Leontief Model (see Economic Models and Problems, Input-Output Systems) 
LIFO vs FIFO (see Inventory and Supply, Qualitative Aspects) 
Linear Equations (see Miscellaneous Mathematics) 
Linear Inequalities (see Miscellaneous Mathematics) 
Linear Programming (see Programming, Linear) 
Logistics, Concepts and Problems 
"“Logistics—What Is It?'' - H. E. Eccles, Vol. 1, No. 1, March 1954, pp. 5-15. 
"Trends in Logistics" - L. R. St. John, Vol. 1, No. 3, September 1954, pp. 182-190. 
"Logistics—The Word" - J..D. Hayes, Vol. 1, No. 3, September 1954, pp. 200-202. 
"Science and Logistics" - R. E. McShane, Vol. 2, No. 1 & 2, March & June 1955, pp. 1-7. 
"Some Command Problems and Decisions" - H. E. Eccles, Vol. 2, No. 1 & 2, March & 
June 1955, pp. 9-15. 
"Note on the Formulation of the Theory of Logistics" - O. Morgenstern, Vol. 2, No. 3, 
September 1955, pp. 129-136. 
"Some General Observations and Experiences in Logistics" - R. B. Carney, Vol. 3, No. 1 
& 2, March & June 1956, pp. 1-9. 
"Concept of a Logistics System" - R. D. Campbell, F. D. Dorey, and R. E. Murphy, 
Vol. 4, No. 2, June 1957, pp. 101-116. 
"Multinational Logistics in the Nuclear Age" - C. L. Henn, Jr., Vol. 4, No. 2, June 1957, 
pp. 117-129. 
“Operations Research and Logistics" - M. M. Flood, Vol. 5, No. 4, December 1958, 
pp. 323-335. 
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Logistics, Historical 
"A New Field for Logistics Research" - B. Anderson, Vol. 1, No. 2, June 1954, pp. 79-81. 
"The Rommel Papers—A Commentary" - H. E. Eccles, Vol. 1, No. 2, June 1954, 
pp. 103-108. 
"The Interdependency of Logistics and Strategic Planning" - S. E. Ruehlow, Vol. 1, No. 4, 
pp. 237-257. 
"Mobile Logistics Support in the ‘Passage to Freedom' Operation" - R. Meade, Jr., and’ 
C. A. Fischer, Vol. 1, No. 4, December 1954, pp. 258-264. 
"Some Command Problems and Decisions" - H. E. Eccles, Vol. 2, No. 1 & 2, March & 
June 1955, pp. 9-15. 
"Some General Observations and Experiences in Logistics'' - R. B. Carney, Vol. 3, No. 1 
& 2, March & June 1956, pp. 1-9. 
Logistics, Planning 
"Some Thoughts on Logistics Planning Factors" - R. E. Williams, Jr., Vol. 1, No. 3, 
September 1954, pp. 173-181. 
"The Interdependence of Logistic and Strategic Planning'"'- S. E. Ruehlow, Vol. 1, No. 4, 
December 1954, pp. 237-257. 
"Consistency Problems in the Military Supply System" - O. Morgenstern, Vol. 1, No. 4, 
December 1954, pp. 265-281. 
"Logistics Without Storage" - I. N. Curtis, Vol. 2, No. 3, September 1955, pp. 125-128. 
"Risk and Hazard in Logistics Planning" - H. A. Sachaklian, Vol. 2, No. 4, December 1955, 
pp. 217-224. 
"Comments on the Morgenstern Model" - J. C. Busby, Vol. 2, No. 4, December 1955, 
pp. 225-236. 
"Improvise or Plan?" - A. H. Wilson and W. R. Finn, Vol. 4, No. 4, December 1957, 
pp. 263-267. 
"Relationships Between Weapons and Logistics Expenditures" - M. A. Geisler, Vol. 4, 
No. 4, December 1957, pp. 335-346. 
"Design Change Impacts on Airframe Parts Inventories" - J. W. Petersen and W. A. 
Steger, Vol. 5, No. 3, September 1958, pp. 241-255. 
Logistics, Survey of Research 
"Some Recent Developments in Transportation Research" - R. R. Crane, Vol. 4, No. 3, 
September 1957, pp. 173-181. 
"Logistics Research Programs of the U.S. Army, the U.S. Air Force and the U. S. 
Navy" - Vol. 5, No. 3, September 1958: 
"Transcript of Briefing on Army Logistics Research" - R. G. Belon, 
pp. 221-223. 
"Briefing on the Logistics Research Program of the Air Force" - R. H. 
Brown, pp. 223-225. 
"Briefing on the Logistics Research Program of the Navy" - C. Stein, Jr., 
pp. 225-229. 
Low Demand (see Inventory and Supply, Demand and Requirements Analysis) 


Machine Scheduling (see Scheduling) 
Maintenance Repair and Salvage 
"The Dynamics of Overhaul and Replenishment Systems for Large Equipments" - 
J. Fennell and S. Oshiro, Vol. 3, No. 1 & 2, March & June 1956, pp. 19-43. 
"Cycling" - I. R. Savage, Vol. 3, No. 3, September 1956, pp. 163-175. 
"On the Theory of Repiacement of Machinery with a Random Failure Time" - G. H. 
Weiss, Vol. 3, No. 4, December 1956, pp. 279-293. 
"Cyclic Scheduling and Combinatorial Topology: Assignment and Routing of Motive Power 
to Meet Scheduling and Maintenance Requirements" - Vol. 4, No. 3, September 1957: 
Part I, "A Statement of the Operating Problem of the Frisco Railroad" - 
V. B. Gleaves, pp. 203-206. 
Part II, "Analysis" - A. Charnes and M. H. Miller, pp. 243-252. 
"Some Queuing Problems in Machine Maintenance" - D. S. Stoller, Vol. 5, No. 1, March 
1958, pp. 83-87. 
"A Transportation and Production Model" - G. W. Evans, II, Vol. 5, No. 2, June 1958, 
pp. 137-154. 
Maintenance and Replacement (see Maintenance Repair and Salvage) 
Management, Inventory (see Inventory and Supply, Management) 












Managerial Economics (see Economic Models and Problems, Managerial Economics) 
Marks=nan vs a Mobile Target (see Games Models) 
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Materiel Improvement Plan (see Economic Models and Problems, Utility) 
Medical 


"Statistical Methods for Determining Requirements of Dental Materials" - J. A. English 
and E. A. Jerome, Vol. 1, No. 3, September 1954, pp. 191-199. 


Miscellaneous Mathematics 


"Computation of Maximal Flows in Networks" - D. R. Fulkerson and G. B. Dantzig, Vol. 2, 
No. 4, December 1955, pp. 277-283. 

"The Basic Theorems of Real Linear Equations, Inequalities, Linear Programming, and 
Game Theory" - D. Gale, Vol. 3, No. 3, September 1956, pp. 193-200. 

"Variants of the Hungarian Method for Assignment Problems" - H. Kuhn, Vol. 3, No. 4, 
December 1956, pp. 253-258. 

"On the Graphs and Asymptotic Forms of Finite Boolean Relation Matrices and Stochastic 
Matrices" - D. Rosenblatt, Vol. 4, No. 2, June 1957, pp. 151-167. 

"Cycling Scheduling and Combinatorial Topology: Assignment and Routing of Motive 
Power to Meet Scheduling and Maintenance Requirements" - Vol. 4, No. 3, Septem- 
ber 1957: 

Part I, "A Statement of the Operating Problem of the Frisco Railroad" - 
V. B. Gleaves, pp. 203-206. 
Part II, "Analysis" - A. Charnes and M. H. Miller, pp. 243-252. 

"Communication-Transportation Networks" (Abstract) - R. E. Kalaba and M. L. Juncosa, 
Vol. 4, No. 3, September 1957, pp. 221-222. 

"Formulation of Recurrence Equations for Shuttle Process and Assembly Line" - 
R. Bellman, Vol. 4, No. 4, December 1957, pp. 321-334. 

"An Optimal Search Pattern" - J. R. Isbell, Vol. 4, No. 4, December 1957, pp. 357-359. 

"Some Queuing Problems in Machine Maintenance" - D. S. Stoller, Vol. 5, No. 1, March 
1958, pp. 83-87. 

"Some Studies on Shuttle and Assembly-Line Processes" - M. Pollack, Vol. 5, No. 2, 
June 1958, pp. 125-136. 

"Nonlinear Network Flows and Convex Programming over Incidence Matrices" - 
A. Charnes and W. W. Cooper, Vol. 5, No. 3, September 1958, pp. 231-240. 

"A Bottleneck Situation Involving Interdependent Industries" - R. Bellman and S. Dreyfus, 
Vol. 5, No. 4, December 1958, pp. 307-314. 

"Multiple Servers with Limited Waiting Space" - J. R. Jackson, Vol. 5, No. 4, December 
1958, pp. 315-321. 

"Network Analysis of Production Smoothing" - T. C. Hu and W. Prager, Vol. 6, No. 1, 
March 1959, pp. 17-23. 

"A Terminating {terative Algorithm for Solving Certain Games and Related Sets of Linear 
Equations" - R. B. Braithwaite, Vol. 6, No. 1, March 1959, pp. 63-74. 

"Prolegomena to Optimum Discrete Search Procedures" - N. M. Blachman, Vol. 6, No. 4, 
December 1959, pp. 273-281. 

"A General Theory of Measurement—Applications to Utility" - J. Pfanzagl, Vol. 6, No. 4, 
December 1959, pp. 283-294. 

"A Matrix Algorithm for Solutions and r-Bases of a Finite Irreflexive Relation" - 
F. Harary and M. Richardson, Vol. 6, No. 4, December 1959, pp. 307-314. 


Military Planning (see Logistics, Planning) 

Military Worth (see Economic Models and Problems, Utility) 
Mission (see Game Theory; also Organization Theory) 

Monte Carlo Methods (see Simulation) 

Multiple Servers (see Miscellaneous Mathematics) 
Multistage Processes (see Scheduling) 


Networks (see Miscellaneous Mathematics) 
Nonlinear Programming (see Programming, Nonlinear) 


Organization Theory 





"On the Span of Central Direction" - T. M. Whitin, Vol. 1, No. 1, March 1954, pp. 25-35. 
"Note on the Formulation of the Theory of Logistics" - O. Morgenstern, Vol. 2, No. 3, 
September 1955, pp. 129-136. 
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Organization Theory —Continued 
"The Concept of a 'Mission''' - M. Verhulst, Vol. 3, No. 1 & 2, March & June 1956, 
pp. 45-57. 
"On the Graphs and Asymptotic Forms of Finite Boolean Relation Matrices and Stochastic 
Matrices" - D. Rosenblatt, Vol. 4, No. 2, June 1957, pp. 151-167. 
Overhaul (see Maintenance Repair and Salvage) 


Parametric Linear Programming (see Programming, Linear) 
Personnel Assignments (see Scheduling; also Programming, Linear—Assignment Problem) 
Planning, Logistical (see Logistics, Planning) 
Policy, Inventory (see Inventory and Supply, Management; also Inventory and Supply, Stocking 
Policies)Portfolio Selection (see Economic Models and Problems, Managerial Economics) 
Primal-Dual Method (see Programming, Linear) 
Priorities (see Economic Models and Problems, Utility; also Inventory and Supply, 
Management) 
Process Analysis, Economic (see Economic Models and Problems, Process Analysis) 
Procurement (see Inventory and Supply, Procurement) 
Production (see Inventory and Supply, Production and Sales) 
Production Scheduling (see Scheduling) 
Productivity Standard, Engineering (see Inventory and Supply, Management) 
Programming, Dynamic 
"On Some Applications of the Theory of Dynamic Programming to Logistics" - R. Bellman, 
Vol. 1, No. 2, June 1954, pp. 141-153. 
"Decision Making in the Face of Uncertainty—I" - R. Bellman, Vol. 1, No. 3, September 
1954, pp. 230-232. 
"Decision Making in the Face of Uncertainty—II"' - R. Bellman, Vol. 1, No. 4, December 
1954, pp. 327-332. 
"The Structure of Dynamic Programming Models" - S. Karlin, Vol. 2, No. 4, December 
1955, pp. 285-294. 
"Notes on the Theory of Dynamic Programming, Part IV, Maximization over Discrete 
Sets" - R. Bellman, Vol. 3, No. 1 & 2, March & June 1956, pp. 67-70. 
"Mathematical Programming and Employment Scheduling"' - W. Karush and A. Vazsonyi, 
Vol. 4, No. 4, December 1957, pp. 297-320. 
"Loss Resulting from the Use of Incorrect Data in Computing an Optimal Inventory 
Policy" - J. Levy, Vol. 5, No. 1, March 1958, pp. 75-81. 
"A Bottleneck Situation Involving Interdependent Industries" - R. Beilman and S. Dreyfus, 
Vol. 5, No. 4, December 1958, pp. 307-314. 
"Further Notes on the Loss Resulting from the Use of Incorrect Data in Computing an 
Optimal Inventory Policy" - J. Levy, Vol. 6, No. 1, March 1959, pp. 25-31. 
Programming, Linear—General 
"The Dual Method of Solving the Linear Programming Problem" - C. E. Lemke, Vol. 1, 
No. 1, March 1954, pp. 36-47. 
"Linear Programming in Bid Evaluation" - E. D. Stanley, D. P..Honig, and L. Gainen, 
Vol. 1, No. 1, March 1954, pp. 48-54. 
"The Caterer Problem" - W. Jacobs, Vol. 1, No. 2, June 1954, pp. 154-165. 
"Minimizing the Number of Tankers to Meet a Fixed Schedule" - G. R. Dantzig and D. R. 
Fulkerson, Vol. 1, No. 3, September 1954, pp. 217-222. 
"On the Solution of the Caterer Problem" - J. W. Gaddum, A. J. Hoffman, and D. Soko- 
lowsky, Vol. 1, No. 3, September 1954, pp. 223-229. 
"The Computational Algorithm for the Parametric Objective Function" - S. Gass and 
T. Saaty, Vol. 2, No. 1 & 2, March & June 1955, pp. 39-45. 
"The Hungarian Method for the Assignment Problem" - H. W. Kuhn, Vol. 2, No. 1 & 2, 
March & June 1955, pp. 83-97. 
"A Dynamic Transportation Model" - J. E. Kelley, Jr., Vol. 2, No. 3, September 1955, 
pp. 175-180. 
"Cycling in the Dual Simplex Algorithm" - E. M. L. Beale, Vol. 2, No. 4, December 1955, 
pp. 269-275. 
"Computation of Maximal Flows in Networks" - D. R. Fulkerson and G. B. Dantzig, 
Vol. 2, No. 4, December 1955, pp. 277-283. 
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Programming, Linear—General—Continued 
"A Plan to Allocate and Procure Electronic Sets by the Use of Linear Programming Tech- 
niques and Analytical Methods of Assigning Values to Qualitative Factors" - J. W. 
Smith, Vol. 3, No. 3, September 1956, pp. 151-162. 
"The Basic Theorems of Real Linear Equations, Inequalities, Linear Programming, and 
Game Theory" - D. Gale, Vol. 3, No. 3, September 1956, pp. 193-200. 
"A Bound on the Use of Inefficient Indivisible Units" - M. J. Beckmann, and J. Laderman, 
Vol. 3, No. 4, December 1956, pp. 245-252. 
"Procurement and Allocation of Naval Electronic Equipments" - G. Sazuki, Vol. 4, No. 1, 
March 1957, pp. 1-7. 
"Linear Program Planning of Refinery Operations" - L. K. Cheney, Vol. 4, No. 1, March 
1957, pp. 9-16. 
"Experiences with the Bid Evaluation Problem" (Abstract) - H. Bremer, W. Hall, and 
M. Paulsen, Vol. 4, No. 1, March 1957, p. 27. 
"A Threshold Method for Linear Programming" - J. E. Kelley, Jr., Vol. 4, No. 1, March 
1957, pp. 35-45. 
"A Primal-Dual Algorithm for the Capacitated Hitchcock Problem" - L.’R. Ford and 
D. R. Fulkerson, Vol. 4, No. 1, March 1957, pp. 47-54. 
"Translating the Method of Reduced Matrices to Machines" - B. A. Galler and P. S. 
Dwyer, Vol. 4, No. 1, March 1957, pp. 55-71. 
"Constraint Matrices of Transportation-Type Problems" (Abstract) - I. Heller, Vol. 4, 
No. 1, March 1957, pp. 73-76. 
"On the Assignment and Transportation Problems" (Abstract) - J. Munkres, Vol. 4, No. 1, 
March 1957, pp. 77-78. 
"Loss of Accuracy in Simplex Computations" - W. Jacobs, Vol. 4, No. 1, March 1957, 
pp. 89-94. 
"Some Methods of Computational Attack on Programming Problems Other than the Simplex 
Method" (Abstract) - C. Tompkins, Vol. 4, No. 1, March 1957, pp. 95-96. 
"Editing Large Linear Programming Matrices" - P. M. Thompson, Vol. 4, No. 1, March 
1957, pp. 97-100. 
"An Analysis of Stewardess Requirements and Scheduling for a Major Domestic Airline" - 
J. F. McCloskey and F. Hanssmann, Vol. 4, No. 3, September 1957, pp, 183-191: 
Annex A, .""The Assignment Problem Technique" - L. Friedman and A. J. 
Yaspan, pp. 193-197. 
Annex B, "Calculation of the School Table" - L. Friedman, pp. 199-202. 
"Mathematical Programming and Evaluation of Freight Systems" - Vol. 4, No. 3, 
September 1957: 
Part I, "Applications" - H. R. Soyster, pp. 237-242. 
Part II, "Analysis'' - A. Charnes and M. H. Miller, pp. 243-252. 
"Numerical Solution of the Generalized Transportation Problem" - W. Prager, Vol. 4, 
No. 3, September 1957, pp. 253-261. 
"On a Certain Problem in Linear Programming" - R. J. Taylor and S. P. Thompson, 
Vol. 5, No. 2, June 1958, pp. 171-187. 
"The Dual Simplex Algorithm for Bounded Variables"' - H. M, Wagner, Vol. 5, No. 3, 
September 1958, pp. 257-261. 
"Operations Research and Logistics" - M. M. Flood, Vol. 5, No. 4, December 1958, 
pp. 323-335. 
"Note on Solving Linear Programs in Integers" - G. B. Dantzig, Vol. 6, No. 1, March 
1959, pp. 75-76. 
"A Difficulty in Linear Programming for Transportation Problems" - M. Luther and 
J. E. Walsh, Vol. 6, No. 1, March 1959, pp. 77-80. 
"Maximum Number of Iterations in the Transportation Problem" - M. A. Simonnard and 
G. F. Hadley, Vol. 6, No. 2, June 1959, pp. 125-129. 
"An Integer Linear-Programming Model for Machine Scheduling" - R. M. Wagner, Vol. 6, 
No. 2, June 1959, pp. 131-140. 
"A Simplified Two-Phase Technique for the Simplex Method" - G.. F. Hadley and M. A. 
Simonnard, Vol. 6, No. 3, September 1959, pp. 221-226. 
Programming, Linear—The Transportation Problem 
"The Dual Method of Solving the Linear Programming Problem" - C. E. Lemke, Vol. 1, 
No. 1, March 1954, pp. 36-47. 
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Programming Linear—The Transportation Problem—Continued 

"Minimizing the Number of Tankers to Meet a Fixed Schedule" - G. B. Dantzig and D. R. 
Fulkerson, Vol. 1, No. 3, September 1954, pp. 217-222. 

"A Dynamic Transportation Model" - J. E. Kelley, Jr., Vol. 2, No. 3, September 1955, 
pp. 175-180. 

"A Threshold Method for Linear Programming" - J. E. Kelley, Jr., Vol. 4, No. 1, March 
1957, pp. 35-45. 

"A Primal-Dual Algorithm for the Capacitated Hitchcock Problem" - L. R. Ford, Jr., 
and D. R. Fulkerson, Vol. 4, No. 1, March 1957, pp. 47-54. 

"Translating the Method of Reduced Matrices to Machines" - B. A. Galler and P. S. 
Dwyer, Vol. 4, No. 1, March 1957, pp. 55-71. 

"Constraint Matrices of Transportation-Type Problems" (Abstract) - I. Heller, Vol. 4, 
No. 1, March 1957, pp. 73-76. 

"On the Assignment and Transportation Problems" (Abstract) - J. Munkres, Vol. 4, No. 1, 
March 1957, pp. 77-78. 

"Mathematical Programming and Evaluation of Freight Shipment Systems" - Vol. 4, No. 3, 
September 1957: 

Part I, "Applications" - H. R. Soyster, pp. 237-242. 
Part II, "Analysis" - A. Charnes and M. H. Miller, pp. 243-252. 

"Numerical Solution of the Generalized Transportation Problem" - W. Prager, Vol. 4, 
No. 3, September 1957, pp. 253-261. 

"A Transportation and Production Model" - G. W. Evans, Il, Vol. 5, No. 2, June 1958, 
pp. 137-154. 

"An Algorithm for Solving the Transportation Problem when the Shipping Cost over Each 
Route is Convex" - E. M. L. Beale, Vol. 6, No. 1, March 1959, pp. 43-56. 

"A Difficulty in Linear Programming for Transportation Problems" - M. Luther and 
J. E. Walsh, Vol. 6, No. 1, March 1959, pp. 77-80. 

"Maximum Number of Iterations in the Transportation Problem" - M. A. Simonnard and 
G. F. Hadley, Vol. 6, No. 2, June 1959, pp. 125-129. 

Programming, Linear—The Assignment Problem 

"The Hungarian Method for the Assignment Problem" - H. W. Kuhn, Vol. 2, No. 1 & 2, 
March & June 1955, pp. 83-97. 

"Variants of the Hungarian Method for Assignment Problems" - H. W. Kuhn, Vol. 3, 

No. 4, December 1956, pp. 253-258. 

"On the Assignment and Transportation Problems" (Abstract) - J. Munkres, Vol. 4, No. 1, 
March 1957, pp. 77-78. 

"An Analysis of Stewardess Requirements and Scheduling for a Major Domestic Airline" - 
J. F. McCloskey and F. Hanssmann, Vol. 4, No. 3, September 1957, pp. 183-191: 

Annex A, "The Assignment Problem Technique" - L. Friedman and A. J. 
Yaspan, pp. 193-197. 
Annex B, "Calculation of the School Table" - L. Friedman, pp. 199-202. 

"The Coefficients in an Allocation Problem" - R. J. Aumann and J. B. Kruskal, Vol. 5, 
No. 2, June i958, pp. 111-123. 

Programming, Nonlinear 

"Optimal Weapon Systems" - J. Marschak and M. R. Mickey, Vol. 1, No. 2, June 1954, 
pp. 116-140. 

"Minimization of Nonlinear Separable Convex Functionals" - A. Charnes and C. E. 
Lemke, Vol. 1, No. 4, December 1954, pp. 301-312. 

"An Algorithm for Quadratic Programming" - M. Frank and P. Wolfe, Vol. 3, No. 1 & 2, 
March & June 1956, pp. 95-110. 

"The Optimization of a Quadratic Function Subject to Linear Constraints" - H. Markowitz, 
Vol. 3, No. 1 & 2, March & June 1956, pp. 111-133. 

"A Quadratic Programming Procedure" - C. Hildreth, Vol. 4, No. 1, March 1957, pp. 79- 
85; also see Errata Sheet, Vol. 4, No. 4, December 1957, p. 361. 

"Computing Procedures for Portfolio Selection" (Abstract) - H. Markowitz, Vol. 4, No..1, 
March 1957, pp. 87-88. 

"The Allocation of Development Funds: An Analytic Approach" - S. C. Daubin, Vol. 5, 
No. 3, September 1958, pp. 263-276; also see Errata Sheet, Vol. 6, No. 2, pp. 
185-186. 

"Nonlinear Network Flows and Convex Programming over Incidence Matrices" - A. Charnes 
and W. W. Cooper, Vol. 5, No. 3, September 1958, pp. 231-240. 
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Programming Nonlinear—Continued 
"Redistribution of Total Stock over Several User Locations" - S. G. Allen, Vol. 5, No. 4, 
December 1958, pp. 337-345. 
"A Convex Programming Procedure" - D. A. D'Esopo, Vol. 6, No. 1, March 1959, 
pp. 33-42. 
"An Algorithm for Solving the Transportation Problem When the Shipping Cost over Each 
Route is Convex" - E. M. L. Beale, Vol. 6, No. 1, March 1959, pp. 43-56. 
"On Quadratic Programming" - E. M. L. Beale, Vol. 6, No. 3, September 1959, pp. 
227-243. 
"A Theorem in Convex Programming" - W. Karuch, Vol. 6, No. 3, September 1959, 
pp. 245-260. 
Pursuit Games (see Games Models) 






Quadratic Programming (see Programming Nonlinear) 
Quality Control (see Reliability and Quality Control) 
Queuing (see Miscellaneous Mathematics) 


Record Storage (see Data Processing) 
Redistribution (see Inventory and Supply, Distribution) 
Refinery Operations (see Economic Models and Problems, Process Analysis) 
Reliability and Quality Control 
"A Mathematical Evaluation of a Work Sampling Technique" - H. Davis, Vol. 2, No. 1 & 2, 
March & June 1955, pp. 111-117. 
"On the Theory of Replacement of Machinery with a Random Failure Time" - G. H. Weiss, 
Vol. 3, No. 4, December 1956, pp. 279-293. 
"The New Military Standard 414 for Acceptance Inspection by Variables" - H. R. Rosen--: 
blatt, Vol. 6, No. 2, June 1959, pp. 173-182. 
Renewal Theory (see Inventory and Supply, General) 
Repair Parts (see Maintenance Repair and Salvage) 
Replacement (see Maintenance Repair and Salvage) 
Research, Surveys of Logistics (see Logistics, Surveys of Research) 











































Salvage (see Maintenance Replacement and Salvage) 
Sampling Methods (see Statistical Decision Theory) 
Scheduling 
"Optimal Two- and Three-Stage Production Schedules with Setup Times Included" - S. M. 
Johnson, Vol. 1, No. 1, March 1954, pp. 61-68. 
"Minimizing the Number of Tankers to Meet a Fixed Schedule" - G. B. Dantzig and D. R. 
Fulkerson, Vol. 1, No. 3, September 1954, pp. 217-222. 
"The Dynamics of Overhaul and Replacement Systems for Large Equipments" - J. Fennell 
and S. Oshiro, Vol. 3, No. 1 & 2, March & June 1956, pp. 19-43. 
"Various Optimizers for Single-Stage Production" - W. E. Smith, Vol. 3, No. 1 & 2, March 
& June 1956, pp: 59-66. 
"Cycling" - I. R. Savage, Vol. 3, No. 3, September 1956, pp. 163-175. 
"An Extension of Johnson's Results on Job Lot Scheduling" - J. R. Jackson, Vol. 3, No. 3, 
September 1956, pp. 201-203. 
"On the Distribution of Manhours to Meet Scheduled Requirements" - S. I. Gass, Vol. 4, 
No. 1, March 1957, pp.-17-25. 
"An Algorithm for the Minimum Number of Transport Units to Maintain a Fixed 
Schedule" - T, E. Bartlett, Vol. 4, No. 2, June 4957, pp. 139-149. 
"An Analysis of Stewardess Requirements and Scheduling for a Major Domestic Airline" - 
J. F. McCloskey and F. Hanssmann, Vol. 4, No. 3, September 1957, pp. 183-191: 
Annex A, ''The Assignment Problem Technique" - L. Friedman and A. J. 
Yaspan, pp. 193-197. 
Annex B, "Calculation of the School Tables" - L. Friedman, pp. 199-202. 
"Cyclic Scheduling and Combinatorial Topology: Assignment and Routing of Motive 
Power to Meet Scheduling and Maintenance Requirements" - Vol. 4, No. 3, 
September 1957: 
Part I, "A Statement of the Operating Problem of the Frisco Railroad" - 
V. B. Gleaves, pp. 203-206. 
Part II, "Analysis" - A. Charnes and M. H. Miller, pp. 243-252. 
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Scheduling —Continued 
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"Communication-Transportation Networks" (Abstract) - R. E. Kalaba, and H. L. Juncosa, 
Vol. 4, No. 3, September 1957, pp. 221-222. 
"Mathematical Programming and Evaluation of Freight Shipment Systems" - Vol. 4, 
No. 3, September 1957: 
Part I, "Applications" - H. R. Soyster, pp. 237-242. 
Part II, "Analysis" - A. Charnes and M. H. Miller, pp. 243-252. 
"Some Production Planning Problems" - M. Klein, Vol. 4, No. 4, December 1957, 
pp. 269-286. 
"Simulation Research on Job Shop Production" - J. R. Jackson, Vol. 4, No. 4, December 
1957, pp. 287-295. 
"Mathematical Programming and Employment Scheduling"' - W. Karush and A. Vazsonyi, 
Vol. 4, No. 4, December 1957, pp. 297-320. 
"Formulation of Recurrence Equations for Shuttle Process and Assembly Line" - R. Bell- 
man, Vol. 4, No. 4, December 1957, pp. 321-334. 
"Some Studies on Shuttle and Assembly-Line Processes" - M. Pollack, Vol. 5, No. 2, 
June 1958, pp. 125-136. 
"Operations Research and Logistics" - M. M. Flood, Vol. 5, No. 4, December 1958, 
pp. 323-335. 
"Network Analysis of Production Smoothing" - T. C. Hu and W. Prager, Vol. 6, No. 1, 
March 1959, pp. 17-23. 
"An Integer Linear-Programming Model for Machine Scheduling" - H. M. Wagner, Vol. 6, 
No. 2, June 1959, pp. 131-140. 
Search Theory (see Miscellaneous Mathematics) 
Seasonal Effects (see Inventory and Supply, Production and Sales) 
Sequencing (see Scheduling) 
Shuttle Processes (see Cargo Handling -and Scheduling) 
Simplex Method (see Programming, Linear) 
Simulation 
"Analysis and Monte Carlo Simulation of Cargo Handling" - R. R. O'Neill, Vol. 4, No. 3, 
September 1957, pp. 223-236. 
"Simulation Research on Job Shop Production" - J. R. Jackson, Vol. 4, No. 4, December 
1957, pp. 287-295. 
Statistics (see Reliability and Quality Control; also Statistical Decision Theory) 
Statistical Decision Theory 
"Sequential Tests of Hypotheses About the Mean Occurrence Time of a Continuous Param- 
eter Poisson Process" - J. Kiefer and J. Wolfowitz, Vol. 3, No. 3, September 1956, 
pp. 205-219. 
"The New Military Standard 414 for Acceptance Inspection by Variables"' - H. R. Rosen- 
blatt, Vol. 6, No. 2, June 1959, pp. 173-182. 
"A Method for Determining Supply Quantity for the Case of Poisson Distribution of 
Demand" - M. R. Mickey, Vol. 6, No. 4, December 1959. 
Stock Levels (see Inventory and Supply, Stocking Policies) 
Stocking Policies (see Inventory and Supply, Stocking Policies) 
Stockpiling (see Inventory and Supply, Stocking Policies) 
Storage (see Records Storage; also Inventory and Supply Warehousing; also Logistics, 
Planning ) 
Strategic Planning (see Logistics, Planning) 
Supply Management (see Inventory and Supply, Management) 
Survival, Games of (see Games Models) 
Systems Analysis 
"Optimal Weapon Systems" - J. Marschak and M. R. Mickey, Vol. 1, No. 2, June 1954, 
pp. 116-140. 
"Trends in Logistics" - L. R. St. John, Vol. 1, No. 3, September 1954, pp. 182-190. 
"System Analysis and/or Common Sense" - Y. Rollof, Vol. 3, No. 1 & 2, March & June 
1956, pp. 11-18. 
"A Formulation of the Decision Problem for a Class of Systems" - J. K. Hale and R. L.. 
Reed, Vol. 3, No. 4, December 1956, pp. 259-277. 
"Relationships Between Weapons and Logistics Expenditures" - M. A. Geisler, Vol. 4, 
No. 4, December 1957, pp. 335-346. 
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Systems Analysis—Continued 
"An Application of Game Theory to Special Weapons Evaluation" - J. K. Hale and H. H. 
Wicke, Vol. 4, No. 4, December 1957, pp. 347-356. 
"Linear Model for Evaluating Complex Systems" - R. M. Thrall, C. H. Coombs, and 
W. Caldwell, Vol. 5, No. 4, December 1958, pp. 347-361. 
Systems, Logistics (see Logistics, Concepts and Problems; also Logistics, Planning) 


Theory of the Firm (see Economic Models and Problems, General) 
Tests of Hypotheses (see Statistical Decision Theory) 
Threshold Method (see Programming, Linear) 
Transportation 
"The Cargo Handling System" - R. R. O'Neill and J. K. Weinstock, Vol. 1, No. 4, 
December 1954, pp. 282-288. 
"A Mathematical Evaluation of a Work Sampling Technique" - H. Davis, Vol. 2, No. 1 & 2, 
March & June 1955, pp. 111-117. 
"A Bound on the Use of Inefficient Indivisible Units'' - M. J. Beckman and J. Laderman, 
Vol. 3, No. 4, December 1956, pp. 245-252. 
"An Algorithm for the Minimum Number of Transport Units to Maintain a Fixed 
Schedule" - T. E. Bartlett, Vol. 4, No. 2, June 1957, pp. 139-149. 
"Communication-Transportation Networks" (Abstract) - R. E. Kalaba and M. L. Juncosa, 
Vol. 4, No. 3, September 1957, pp. 221-222. 
"Analysis and Monte Carlo Simulation of Cargo Handling" - R. R. O'Neill, Vol. 4, No. 3, 
September 1957, pp. 223-236. : 
"Mathematical Programming and Evaluation of Freight Shipment Systems" - Vol. 4, No. 3 
September 1957: 
Part I, "Applications" - H. R. Soyster, pp. 237-242. 
Part II, "Analysis" - A. Charnes and M. H. Miller, pp. 243-252. 
"Formulation of Recurrence Equations for Shuttle Process and Assembly Line" - R. Bell- 
man, Vol. 4, No. 4, December 1957, pp. 321-334. 
"Some Studies on Shuttle and Assembly-Line Processes" - M. Pollack, Vol. 5, No. 2, 
June 1958, pp. 125-136. 
"Characteristics of Usage of Supply Items Aboard Naval Ships and the Significance to 
Supply Management" - J. H. Garrett, Jr., Vol. 5, No. 4, December 1958, pp. 287-306. 
"Operations Research and Logistics" - M. M. Flood, Vol. 5, No. 4, December 1958, pp. 
323-335. 


"The Feasibility of a Global Airlift - R. S. Weinberg and J. W. Higgins, Vol. 6, No. 2, 
June 1959, pp. 89-110. 
"Allocation of a Resource to Alternative Probabilistic Demands: Transport-Equipment 
Pool Assignments" - R. Saposnik and V. L. Smith, Vol. 6, No. 3, September 1959, 
pp. 193-207. 
Transportation Problem, The (see Programming, Linear—The Transportation Problem) 
Transportation Scheduling (see Scheduling) 


Usage Data (see Inventory and Supply, Demand and Requirements Analysis) 
Utility (see Economics Models and Problems, Utility) 
Utility Theory (see Miscellaneous Mathematics) 


Warehousing (see Inventory and Supply, Warehousing) 

War Games (see Games Models) 

Waiting Line Problems (see Scheduling; also Miscellaneous Mathematics) 
Work Sampling (see Reliability and Quality Control) 

Weapons Evaluation (see Systems Analysis) 

Weapons Systems (see Systems Analysis) 
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